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APPENDIX  A 


DOA  ACCURACY  ANALYSIS 


One  technique  for  estimation  of  directlon-of-arrival  (DOA)  of 
acoustic  or  electromagnetic  energy  is  the  use  of  sensors  or  antennas 
whose  sensitivity  is  a  function  of  sine  and  cosine  of  DOA  of  the  incident 
energy.  The  signals  received  with  an  omnidirectional  sensor  and  two 
such  sinusoidal  sensors  (for  a  tone  in  noise)  are  given  below: 


e0(t) 

es(t> 

ec(t) 


VB(t)cos  [cat  +  0(t)]  +  nQ(t) 
VB(t)sin9  cos  [cot  +  0(t)]  +  ng(t) 
VB(t)cos  9  cos  [cot  +  0(t)]  +  n£(t) 


(A-l) 


(A-2) 


(A-3) 


where  n  (t) ,  n  (t) ,  and  n  (t)  are  assumed  zero  mean  Gaussian  noise  terms , 
OS  c 

0(t)  is  a  uniformly  distributed  (0  to  2ir)  random  variable,  and  B(t)  is  a 
normalized  Rayleigh  distributed  fading  factor.  The  expected  value  of  B(t) 
is  \J ^  and  the  expected  value  of  the  noise-free  envelope  of  the  output  of 
the  omni  channel  is  V .  The  processing  assumed  for  extraction  of  the 
bearing  information,  9,  involves  the  cross-correlation  of  the  omni  channel 
with  both  the  sine  and  cosine  channels  which  removes  terms  at  w  and  2w  and 
yields: 


A-l 


•os™ 


-V2B2(t)sin0  +  — VB(t)n  .  (t)  + -— -VB(t)n  .(t)sin  G 

11  /2?  81  /2? 


+  Ki(t)  nsi(t)  +  Kq(t)  nsq(t) 


(A-4) 


e  (t)  =  —  V2B2(t)cos0  +  — - — VB(t)n  .  (t)  +  — —  VB(t)n  (t)cos0 

oc  *  sn  ci  su  01 


+  2noi(t)  nci(t)  +J  noq(t)ncq(t) 


(A-4) 


,  The  noise  terms  follow  from  an  expansion  of  the  input  noise  into 
their  in-phase  and  quadrature  components,  that  is 

(A-6) 
(A- 7) 


nQ(t)  =  noi(t)cosi  t  +  0)  +  n  (t)  sin(  t  +  0) 


n  (t)  =  n  .  (t)cos(  t  +  0)  +  n  (t)  sin(  t  +  0) 

8  SI  sq 


n^(t)  *  nci(t)cos(  t  +  0)  +  n _ (t)  sin(  t  +  0) 


cq 


(A- 8) 


If  we  assume  that  the  noise  field  is  Isotropic,  it  is  easy  to  show  that 


the  terms  n  (t) ,  n  (t) ,  and  n  (t)  are  mutually  uncorrelated  as  follows:  The 
OS  c 


noise  terms  may  be  defined  as  the  integrals  of  the  noise  field  over  the  full 
360°  of  azimuth  (we  are  assuming  here  that  the  noise  field  is  adequately 
represented  by  one  dimension) : 


2ir 


no(t)  =/  vo(0’t) 


d0 


(A-9) 


The  cross  product  n  (t)n  (t)  can  then  be  written  in  terms  of  v  (0,t): 

OS  o 


2n  2  it 


no(t)ng(t)  -/■/  v(0^ ,t)  vq  (02,t)  sin  dQ^  d0j 


(A-10) 


0  0 


!  1 


I 

J 


A-2 


The  expected  value  of  this  expression  can  be  written  as: 


E[0o<t)ns(t)]  E[Vo^6l,t^Vo^02,t^]  sln6id0i 

0  0 


d0„ 


(A-ll) 


In  general,  noise  terms  originating  from  different  azimuthal  bearings  will 
be  uncorrelated  (discounting,  for  example,  focused  back-scatter  from  a 
rough  bottom).  With  the  noise  assumed  stationary,  the  expected  value  can 

be  written  in  terms  of  a  delta  function: 

2ir  2ir 

E;£no(t)ns(t)J  =  J  J  N(61)6(01-02)  Sin91d01d02 
0  0 


2ir 

= J  N(0)  sin9d0  (A-12) 

0 


For  the  isotropic  noise  field  assumed  here,  the  azimuthal  noise  power  density 
is  independent  of  0  and  can  be  brought  outside  the  integral  as  a  constant: 

2  IT 

E|no(t)ng(t)J  =  N  J  sin  0  d0  =  0 
0 

Similar  arguments  show  that  for  an  isotropic  noise  field  the  omni,  sine, 
and  cosine  channel  noise  terms  are  all  mutually  uncorrelated  (and,  there¬ 
fore,  statistically  independent  under  a  Gaussian  noise  assumption). 

With  these  results,  the  expected  values  of  e  (t)  and  e  (t)  are  readily 

OS  oc 

determined : 


[>“>] 

-  (l/n)V2  sin  0  E  [] 

B2(t)] 

(A-13) 

M>] 

=  (l/jt)V2  cos  0  E  £] 

b2U)] 

(A-14) 

A-3 


Under  the  Rayleigh  fading  assumption,  the  expected  value  of 
equals  2  and  these  relations  can  be  rewritten: 


E[eos(t)]  *=  (2/*)V  sin  0 
■  [•ocCO]  "  (2/«)V2cos  0 


(A-15) 

(A-16) 


Also  under  the  Rayleigh  fading  assumption,  the  expected  value  of  equals 
8  and  the  variance  of  these  "signal"  terms  are: 


e{[(1/it)V2B2  sin0]2  -  [(2/ir)V2sin  0]  2 J 


(A -17) 


a  =  (2/ir)v  | sin  0  | 


(A-18) 


The  similar  result  for  the  cosine  channel  is: 


oc  ■  (2/ir)  y  Jcos  0 1 

With  these  results,  it  is  convenient  to  rewrite  e  (t)  and  e  (t)  in  terms 

os  oc 

of  a  random  variable  e(t)  with  zero  mean  and  unit  variance: 


e  (t)  =  (2/ir)  V  sin  0  <1  +  €  (t))  +  n  (t) 

OS 


(A- 19) 


e  (t)  =  (2/ir>V  cos'  6  (1  +  €  (t))  +  n  (t) 
oc  yt* 


(A-20) 


where  e  (t)  a  B  ^  ^ 

r 

nos(t)  -  VB(t)[n  (t)  +  n  (t)sin  «]//£  +  [n  (t)n  (t)  +  n  (t)n  (t)]  fl 


(A-21 


noc(t)  "  VB(t)  fncl (t)  +  nQl(t)  cos  0]  /yfcx  ■+  lnoi(t)nci(t)  +  noq(t)ncq(t)V2 


A-4 


(A-22) 


The  noise  terns  n  (t)  and  n  (t)  clearly  have  zero  means,  since  n  , 
os  •  oc  o 

and  n  are  uncorrelated.  As  is  well  known,  the  quadrature  components 
8 

are  also  uncorrelated,  if  a  symmetric  noise  spectrum  is  assumed.  The  var< 
iance  can  be  determined  by  finding  the  expected  value  of  their  squares: 


l 

I 


<4  •  e[(”oS)2] 

+(4oi  Mv  ] +  4VM4  ])/“ 


-  E[(v>2] 

■  v2e[b2](E["L  ]  +  E[”ol  Jc08?e)/2' 

♦  (E[noi  JE[nci  ] +  4’J42J)/‘  <a 

It  can  be  shown  that  the  variances  of  the  in-phase  and 
quadrature  noise  amplitude  components  equal  the  variance  of  the  entire 
noise  process,  that  is 

4o2)-  4i2l  ■  E[no,2]  -%  <A 

482)-M -[>]  -  »2  <a 

E(“c2)  ■  E["cl2]  '  E[ncq2]  ■  °C  <A 

These  values  can  be  determined  by  extending  the  results  of  Equation  A-12. 
2 

The  variance  oq  results  by  replacing  ng(t)  with  n  (t).  This  eliminates 

2 

the  sine  function  in  the  integral  and  yields  oq  ■*  2itN.  Replacing  nQ(t) 

2 

with  ng(t)  in  Equation  A-12  would  yield  og.  The  integral  would  contain 


2 

o 

oc 


A-5 


2  2  2 
a  sine  function  and  yield  a  ■  2nN/2.  Similarly  a  can  be  shown  equal 

B  C 

2 

to  0  .  Therefore,  for  an  isotropic  noise  field: 

8 

o2  «=  o2  =  2rtN  (A-27a) 

a2  «  o2/2  (A-27b) 

s 

a2  =  o2/2  (A-27c) 

c 

These  expressions  can  be  used  to  evaluate  Equations  A-23  and 

A-24: 

°os  "  2y2  [(°2/2)  +  °2  sin2  ®]/2ir  +  [(o^2/2)  +  (oZo2/2)J/4 

=  V2o2(1  +  2  sin2  0)/2tt  +  a4/4  (A-28) 

o2  -  V2a2(l  +  2cos2e)/2ir  +  0V4  (A-29) 

oc  •  ' 


The  numbers  s(M)  and  c(M)  result  from  averaging  M  independent  samples  of 

the  processes  e  (t)  and  e  (t),  respectively: 
os  oc 


S(M)  *  (2/^)^  sin  0(1  +  eM)  +  ngM 

(A-30) 

C(M)  =  (2/n)V2  sin  0(1  +  eM>  +  n^ 

(A-31) 

where 

E[£tl]  '  E["sm]’  E[”ch]‘  0 

(A-32a,b ,c) 

e[4]  -  1/M 

(A-32d) 

E[”L]  '  °liK 

(A-32e) 

-  °L/m 

(A-32f ) 

A-6 


To  estimate  the  DOA,  9,  the  ratio  S(M)/C(M)  is  formed* 


S(M)/C(M)  = 


^(2/7r)v2sine(l  +  eM)  +  ngMJ^J(2/Tt)v2cos9(l  +  eM)  +  ncMJ 


(A-33) 


In  the  noise-free  case  this  expression  reduces  to  sin9/cos9  or 
tan9  which  can  be  used  directly  to  estimate  signal  bearing  angle.  In  the 
presence  of  noise,  the  use  of  this  tangent  relationship  will  lead  to 
errors  in  estimation  of  the  bearing  angle  9. 

A  number  of  authors  have  studied  the  errors  associated  with  a 
bearing  estimator  such  as : 


&  =  tan"lj^s(M)/C(M)J 


(A -34) 


In  particular,  Richter  (Ref.  A-l)  studied  the  following  estimator: 


G- 


tan 


1 


Asin9  + 


Acos9  + 


ny(H) 


nx(M) 


(A-35) 


His  results  can  be  applied  to  the  estimator  of  Equation  A-34 
UBing  a  redefinition  of  terms  appropriate  to  our  problem: 

A  -  2V/tt  (A-36) 

VH>  -  "SM/(1  +  'M  )  <A-37) 

nx(M)  "  +  *M)  <A-38) 

In  writing  Equations  A-37  and  A-38  we  recognize  the  fact  that, 
whereas  the  terms  (1  +  e  )  are  associated  with  the  signal  components  in 
forming  S(M)  and  C(M) ,  the  effect  of  this  term  is  equivalently  an  inverse 
scaling  of  the  noise  terms  ng^  and  n^ 


A -7 


In  order  to  apply  Richter's  results,  the  terms  n  (M)  and  n  (M) 

x  y 

should  be  independent  Gaussian  variables  We  notice  that,  since  n  and 

n^  are  uncorrelated,  nx(M)  and  ny(M)  are  also  uncorrelated.  For  large  M, 

the  terms  n  n  and  f.,  may  also  be  assumed  to  be  normally  distributed.  The 
sM  CM  M 

ratios  used  to  define  n  (M)  and  n  (M) »  however,  are  not  normally  distribu- 

y  x 

ted.  For  large  M,  the  variance  of  f  becomes  quite  small  and  in  the  limit, 

as  M  — •>  oo*  goes  to  zero.  In  the  limit,  therefore,  both  n  (M)  and  n  (M) 

y  x 

approach  a  normal  distribution.  In  the  analysis  that  follows  we  will 
assume  M  is  large  enough  for  the  discripancy  to  be  negligible  and  will 

assume  Richter's  results  apply  exactly. 

The  pertinent  portion  of  Richter's  paper  is  presented  in 

2  2 

Appendix  B.  In  order  to  use  Richter's  results,  two  terms,  S  and  a  , 


must  be  defined: 


‘M-i +  %) 


2  2.  2 

i  -  a  /a 

y  x 


where 


a1  =  eL2]  =  eL2  1  eU  +  e  f 

x  xj  L  cMJ  \  «/ 

o2  =  E  n2  =  E  n2M  E  1/ (l  +  e  )2 
y  yj  sMJ  [_  V  M/ 

For  small  eM  (i.e.  large  M)  the  term  E  Hi+e«)2] 
can  be  estimated  by  expanding  1/  ^  1  x 

E[1/(1  +  €m)2]  "  E[X  '  2CM  +  3CM  +  ° (C 3)1 


(A -39) 


(A -40) 


(A-41) 


(A-42) 


(A-43) 


2  2  2 

Substituting  cr  ,  <r  (using  Equations  A-28,  A-29  and  A-43)  and  K  ,  into 
x  y 

Equations  A-39  and  A-40  yields 


(mTs)  (4)’  {*¥  *  f) 


(A-44) 


2  m  V2o2(l  +  2  sin2  9)  /2tt  +  g4/ 4 
V^o2(l  +  2  cos2  9)/2ir  +  a^/4 


(A-45) 


Using  Eq.  A-l,  the  signal-to-noise-ratio  in  the  omni  channel  can 
be  defined: 

Y  =  (2V2/tt)/o2  (A-46) 

This  expression  can  be  used  to  simplify  Equations  A-44  and  A-45: 

.2 


(A-47 


[rd  +  2  sin2  9)  +  l]/[v(l  +  2  cos2  9)  +  l]  (A-48) 


In  general,  M  (the  required  number  of  samples  averaged)  is  the  parameter  which 
determines  system  integration  time.  This  parameter  can  be  determined  from 
Equation  A-47: 


M  -  f  [l  +  \A  +  ^J 

H,f[l+(1  +  |)] 


(A-49) 


M  ~  r  +  3 


(A-50) 


where 


r  =  S2(l  +  2rj/^272 


[  1 


For  large  M,  Ms  r  and  the  condition  necessary  for  the  square  root 
approximation  becomes 


12  <  i 

M  3 


M  >  36 


For  M  =  36,  Equation  A-49  is  within  10%  of  the  exact  expression 
Equation  A-48,  with  the  error  decreasing  to  0  for  large  M.  Also  for 
M  >  36,  Equation  A-43  is  a  very  tight  estimate  of  the  expected  value  of 


Equation  A-49  and  the  curves  of  Appendix  B  were  used  to  generate 
the  curves  of  Figure  A-l,  which  give  rms  bearing  error  as  a  function  of  the 
expected  input  omni-directional  SNR  and  the  number  of  samples  M  for  both  the 

maximum  error  bearing  directions  (45°,  135°,  225°,  315°)  and  the'minimum 
error  bearing  direction  (0°,  90°,  180°,  270°).  In  addition,  since  Richter's 
curves  were  not  readable  for  the  case  of  very  small  error,  a  small  noise 
approximation  was  derived  using  Equations  A-33,  A-34,  and  A-35: 


M  =  [ 


3  -  cos  40 


(A-51) 


where  is  the  desired  rms  bearing  accuracy.  This  expression  was  used  to 
generate  the  1°  rms  curves  of  Figure  A-l.  As  an  example  of  the  derivation 
of  these  curves  and  the  use  of  the  curves  of  Appendix  B,  the  following  set 
of  parameters  are  used  to  determine  M: 

cr 0  =  10°  rms,  (cr^  =  0.03  rad^) 

y  -  -5  dB,  (y  =  0.316) 

2  .  _o 

°0  “  45  >  0 


A- 10 


I 


The  first  parameter  to  be  calculated  is  a.  Using  Equation  A-45, 

it  is  easily  shown  that  at  45°  a  =  l  and  at  0°  a  =  0.82.  With  these  values 
2 

of  Oq  and  a.  Figures  3-2a  (page  B-8)  and  4-2a  (page  B-9)  can  be  used  to 
.  2 

determine  S  .  The  resulting  values  are: 

S2  -  12.5  dB  (17.78),  0  =  45° 

S2  =  11.5  dB  (14.12),  0=0° 

Using  these  values  in  Equation  A-46  yields  the. desired  values  of  M: 

M  =  147  samples,  0  =  45° 

M  =  117  samples,  0=0° 

Finally,  Figure  4-3  (page  B-13)  and  Table  4-1  (page  B-6)  can  be  used  to 
estimate  the  bias  error  resulting  from  S  at  0  =  45°  and  0°.  These  figures 
show  the  resulting  bias  to  be  negligible  at  0°  and  45°  (+njt/2)  . 

For  this  same  set  of  parameters,  the  approximate  expression.  Equation  A-47, 
yields : 

M  =  136  samples,  0  =  45° 

M  «  110  samples,  0=0° 

The  agreement  between  the  numerically  derived  results  and  the  approximate 
results  at  10°  rms  (better  than  8%)  suggests  that  the  approximate  expression 
is  very  tight  at  1°  rms  and  for  many  cases  of  practical  interest  may  be  used 
without  significant  error  up  to  10°  rms. 


A-12 


it  -» 


APPENDIX  B 


AVERAGING  TO  REDUCE  VARIANCE  IN  BEARING  ESTIMATES 


This  appendix  is  a  reproduction  of  part  of  Section  4  and 
Figure  3-22  of  a  1969  Magnavox  Technical  Report  (MX-TR- 6-7 -200169)  of  the 
same  title  which  was  prepared  by  R.  J.  Richter  and  G.  E.  Gaerte.  The 
report  derives  accuracy  expressions  that  result  from  the  use  of  x(t)  and 
y(t)  to  estimate  bearing.  These  are  defined  by 
x(t)  s  A  cos  9  +  n^(t) 
y(t)s  A  sin  9  +  n^(t) 

where  9  is  the  true  signal  bearing.  The  actual  bearing  estimator  analyzed 


0  =  tan-1(y/x) 

In  the  analysis,  the  noise  terms  are  assumed  to  be  zero  mean,  Gaussian, 

2  2 

ergodic  processes  with  unequal  variances  a  and  o  .  In  the  sections 

x  y 

reproduced  here,  three  equations  from  other  sections  are  referenced,  1-30, 
2-21,  and  2-22.  Equation  1-30  is  identically  Equation  4-1  of  the  enclosed 
section.  Equations  2-21  and  2-22  define  the  mean  and  variance  of  0  and  are 
repeated  here: 


*  Jl 

E [ 0 3  - J  0P0(0)d0 

e 

9  +  it 

var[0]  * J  0  P0(0)d0  -  EZ[0] 


(2-21) 


(2-22) 


where  P0  is  the  probability  density  function  of  the  estimator  0. 


B-l 


There  are  also  two  typographical  errors  in  the  material  reproduced. 
The  exponent  of  the  second  exponential  expression  in  Equation  4-1  should 
be  negative  (i.e.,  "-S  (1  +  a  )...")  and  the  exponential  term  in  Equation  4-16 
should  be  scaled  by  a  factor  of  1/2. 


4.0 


CASE  OF  UNEQUAL  VARIANCES,  ZERO  CORRELATION 


Now  consider  the  case  when  X  and  Y  are  uncorrelated  with  un- 

A 

equal  variances.  Again  we  will  first  consider  the  estimator  6  for 
n  =  1.  The  density  function  for  0  is  given  by  equation  (1-30) 
a 


P*(<f>)  = 


-S2(l  +  a2) (a2  cos20  +  sin20) 
exp  { - } 


2ir(a2  cos 2<p  +  sin2$>)  2a2 

S  /l  +  a2  (a2  cos  ip  cos  0  +  sin  <p  sin  6) 


2  /2ir  (a2  cos24>  +  sin24>) 3/2 
S2(l  +  a2)  sin2 (<j>  -  6) 


exp  {- 


j  [1  +  erf  (X)] 


2(a2cos2<f>  +  sin240 
JY  =  ouX 


<jy2(1  +  a2) 


S  /l  +  a?-  (a2  cos  <p  cos  0  +  sin  p  sin  0) 


(4-1) 

(4-2) 

(4-3) 

(4-4) 


X  = 


a  /2  /a2  cos2$  +  sin2$ 

The  analysis  for  this  case  parallels  that  given  in  the  preceding 
paragraphs  for  the  case  of  correlated  signals.  Consequently,  we  will  be 
brief  in  our  discussion  for  this  case. 

4,1  Symmetry  of  Density  Function 

One  can  show  that 


P*(-<!>;  a,  S2,  -9)  =  P*($;  a,  S2,  0) 


(4-5) 


and 


P^OO0  -  a,  S2,  90°  -  0)  =  P4(90°  +  a,  S2,  90°  +  6)  (4-0) 
This  gives  symmetry  about  0=0°  and  0  =  90°  so  that  we  have 

eC-e)  =  -3(0)  (4-7) 

3(90°  ♦  8)  =  -8(90°  -  9)  •  (4-8) 


and 


Var[<i>|-9]  =  Var[$j0] 

Var(90  «•  4>|90  +  0]  Var[90  -  $|90  -  0] 


(4-9) 

(4-10) 


For  example,  the  variance  of  9  for  6  =  -30“  equals  the  variance  for 
e  =  30°,  and  the  varian.ee  of  9  for  0  =  120°  equals  the  variance  for 


6  =  60°. 


One  can  also  show  that 


P*(*  +  90°;  I,  S2,  9  +  90°)  =  P0(*;  a,  S2,  0) 


This  says  that 


Var[4> 1  — ,  0]  =  Var[$|a,  0  -  90°) 


(4-11) 


(4-12) 


For  example,  the  variance  of  9  when  a  =  2  and  9  =  30°  equals  the  variance 
when  a  =  0.5  and  0  =-60“  (or  0  =  60°). 

The  above  symmetries  imply  that  we  need  only  to  consider  the 
density  function  for  0  0  90°  and  0  <  a  <_  1.  Considerations  for  other 

values  of' the  parameters  can  be  made  using  the  symmetry  properties.  Thus 
we  will  restrict  our  attention  to  values  of  0  in  the  first  quadrant  and 
values  of  a  less  than  one. 

4.2  Special  Cases 

4.2.1  S2  »  1 


For  S2  >>  1  and  1 4>  -  0|  <  5°, 
cos  ($  -'  0)  s  1 
sin  (4>  -  0)  s  <j>  -  0 
erf  (X)  =1 


(4-13) 

(4-14) 

(4-15) 


P*(4>;  S2,  a,  0)  = 


27  fa2  cos20  +  sii 
K  S2 (1  +  a2) 


-(*  -  e)?- 

x  exp  { - } 

a2  cos20  +  sin20 

S2(l  +  a2) 


(4-16) 


In  the  vicinity  of  0,  $  is  approximately  Gaussian  with  mean  G  and 

2  2*2  a 

variance  - — ~s- -•  S1.n . One  would  thus  expect  the  estimator  6  to  be 

S2(l  +  a2) 

unbiased  for  large  SNR.  The  variance  is  dependent  upon  a,  S2,  and  G; 
it  is  periodic  in  9  with  period  it.  The  maximum  variance  now  occurs  at 
9  *  +90°  and  the  minimum  occurs  at  0  =  08,  180°.  Note  that  this  is  a 
45°  shift  from  the  case  of  correlated  signals.  The  variance  lies  within 
the  range  a2/S2(l  +  a2)  to  1/S2 ( 1  +  a2)  for  all  9.  The  width  of  this 
range,  (1  -  a2)/S2(l  +  a2),  increases  with  decreasing  a.  The  high  SNR 
asymptote  is 

a2  cos20  +  sin20 

10  log10  S2  =  -10  log10  a2  +  10  log10  ( - }  (4-17) 

1  +  a2 


4.2.2  S2  =  0 

For  the  case  of  no  signal  (S2  =  0)  the  density  function  becomes 

V 

0,'a,  0)  «  - - - 1 - (4-18) 

2n(a2  cos2$  +  sin2<f>) 

The  density  function  is  periodic  with  period  tt  and  has  extreme  points 
at  9  *  0,  +90°,  180°,  As  a  +  1,  Pj  +  2^  which  is  a  uniform  distribution 
with  mean  9  and  variance  tt2/3.  As  a  ■+  0,  the  density  approaches  impulses 
centered  at  $  =  0  and  <j>  =  180°.  We  would  again  expect  the  estimator  to 
be  biased  when  the  SNR  is  low.  Values  of  the  mean  and  variance  for 
various  values  of  a  are  listed  in  Table  4-1. 

4.3  Numerical  Evaluations 

Equations  (2-21)  and  (2-22)  were  evaluated  numerically  for 
several  combinations  of  parameters.  The  results  of  these  evaluations 
are  shown  in  Figures  4-1  to  4-3.  Figure  4-1  shows  several  density 
functions  for  S2  *»  5  dB  and  a  =  0.5.  Typical  plots  of  SNR  versus 
variance  are  shown  in  Figure  4-2.  The  variance  decreases  with  increasing  SNR 
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For  SNR  greater  than  6  dB  the  variance  is  maximum  for  0  =  90°  and 
minimum  for  0=0°.  Near  S2  =  6  dB,  the  variance  is  nearly  independent 

of  0.  •" 

The  estimator  is  biased  for  all  values  of  0  except  0=0,+  90°, 
180°.  Figure  4-3  is  a  plot  of  the  mean  versus  SNR.  The  bias  increases 
with  decreasing  SNR;  the  limiting  values  are  given  in  Table  4-1.  For 
a  >  0.1  the  bias  is  less  than  7°  and  for  a  >  0.5  the  bias  is  less  than 

5° 

°  •  | 

4.4  Analysis  of  Estimators 

The  analysis  of  the  estimators  given  in  paragraph  3.5  also 
applies  to  this  case  of  unequal  variances.  The  discussion  and  test 
described  in  paragraph  3.6  apply  also.  For  high  SNR  both  methods  give 
equal  performance.  The  high  SNR  asymptotes  and  test  points  Pi  and  P2 
are  included  in  Figure  4-2.  From  these  graphs  we  see  that  there  is  again 
a  range  in  which  the  better  method  depends  upon  n.  For  example,  if 
a  =  0.8,  we  see  that  method  2  is  better  for  S2  <  -3  dB  and  method  1  is 
better  for  S2  >  -3  dB  if  n  >  20.  Again  we  see  that  the  better  method 
depends  upon  0  as  well  as  upon  a  and  S2  when  a  <  0.7.  Now  for  many 
practical  considerations  we  will  only  be  concerned  with  values  of  a  >  0.5 
and  S2  >  -5  dB.  In  these  situations  method  1  would  be  preferred  since 
it  is  the  better  method  for  most  values  of  0.'  Also  the  observations 
that  in  general  02  will  not  be  distributed  on  [u  -  0,  n  +  0]  and  that 

A 

the  bias  of  02  does  not  decrease  with  n  tend  to  make  method  1  the 
preferred  method. 
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TABLE  4-1 

NO  SIGNAL  CASE  (S2  =  0) 


8 

E[*]. 

Bias 

Var [4] 

a 

(Degrees) 

(Degrees) 

(Degrees) 

(Rad2) 

0.8 

0 

0.002 

-  0.002 

3.343 

15 

16.581 

-  1.581 

3.335 

30 

32.681 

-  2.681 

3. 314 

45 

48.013 

-  3.013 

3.286 

60 

62.543 

-  2.543 

3.261 

75 

76.441 

-  1.442 

3.244 

90 

89.999 

0.001 

3.238 

0.5 

0 

0.001 

-  0.001 

3.656 

15 

28.188 

-15.188 

3.539 

30 

49.108 

<-19. 108 

3.327 

45 

63.436 

-18.436 

3.159 

60 

73.898 

-13.898 

3.054 

75 

82.369 

-  7.370 

2.998 

90 

90.000 

0.000 

2.9S1 

0.1 

0 

-  0.002 

0.002 

4.401 

15 

69.535 

-54.535 

3.138 

30 

80.174 

-50.174 

2.788 

45 

84.290 

-39.290 

2.674 

60 

86.697 

-26.697 

2.625 

75 

88.466 

-13.466 

2.603 

90 

90.000 

0.000 

2.597 

APPENDIX  C 


SPLIT-BEAM  TRACKER  ACCURACY  ANALYSIS 


A  theoretical  analysis  of  the  performance  of  a  split  beam  tracker 
is  contained  in  Reference  C-l.  One  result  of  that  analysis  is  an  expression 
for  the  rms  pointing  error  for  a  linear  array  of  M  transducers  with  inde¬ 
pendent  noise  outputs  (e.g.  half  wavelength  spacing  between  elements). 


r  w. 


dM  vf”  cos  6 


/ 


do)  (o 


S2(gj)/N2(6>) 

1  +  M[S  (6j)/N(oj)] 


1 

2 


(C-l) 


where,  d  is  the  element  space 

c  is  the  propagation  velocity 

T  is  the  system  integration  time 

M  is  the  number  of  sensors 

0  is  the  azimuthal  arrival  angle  relative  to  broadside 

S(e>)  is  the  signal  power  spectral  density 

N(&>)  is  the  noise  power  spectal  density 

6j  is  the  input  bandwidth  processed  in  radians 
N 

For  systems  not  operating  at  baseband,  the  integration  can  be 
be  changed  from  the  range  0  to  to  the  range  /  ^) 

ft»0+  where  (oQ  is  the  center  frequency.  We  now  assume  that  the  sign 

to  noise  power  spectral  density  ratio  is  constant  over  the  processed  band¬ 
width.  Denoting  this  signal-to-noise  ratio  by  7,  the  integral  is  readily 
evaluated  yielding 


* 


<mVt  cos  e 


ifl2oj2w  + 

Von  n  / 

12(1  +  M7) 


1 

2 


(C-2) 
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Note  that  this  signal  to  noise  ratio  is  that  defined  for  an  omni-directional 


sensor . 


For  narrow  band  signals  (i.e. 
This  expression  simplifies  to: 


Q*  a 


(C-3) 


where  0  is  assigned  near  boresight  and,  therefore,  cos  0  has  been  approximated 
as  unity. 

This  expression  can  be  rewritten  in  terms  of  the  array  length 
L,  the  wavelength  A,  and  the  bandwidth  F,  which  are  defined  by 


L  ■= 
A  = 

F  = 


(M-l)d 

2jtcAun 


CUN 

~  Hertz 

2it 


We  then  have 


For  M  »  1  and  MY  »  1,  this  expression  can  be  simplified  further. 


or  => 


JL-v/H 

jtL  VTFMY 


(C-5) 


Finally,  the  ratio  X/L  can  be  related  to  the  beamwidth,  B>  normal 
to  the  linear  array  (Reference  C-2  page  25-40) : 


B, .  -  0.88A/L 

lin 

This  yields 


cr  = 

0.88jt VTFMY 


(C-6) 


(C-7) 


Under  the  assumption  that  the  environmental  noise  is  isotropic  and,  there¬ 
fore,  independent  from  sensor  to  sensor  at  a  1/2  wavelength  spacing,  M  can 
be  expressed  in  terms  of  the  broadside  beamwidth  using  Equation  C-6: 

M  ~  M- 1  •  2L/A  =  2 

v  lin 


) 


(C-8) 


Substituting  into  Equation  C-7  and  evaluating  the  numerical  form  results 
in  the  desired  expression  for  the  rms  bearing  error  obtained  with  a  linear 
array: 


°-38(Blin) 
VTF7 


3/2 


(C-9) 


We  now  consider  the  use  of  a  circular  planar  rather  than  a 
linear  array.  An  approximate  expression  for  the  estimation  error  result¬ 
ing  from  the  use  of  a  circular  planar  array  can  be  generated  by  recogniz¬ 
ing  that  in  an  infinite  ocean  of  finite  depth  the  vertical  beam  pattern 
provides  little  gain  over  an  omni-directional  sensor  since  both  the  omni 
and  vertically  directional  patterns  receive  energy  from  the  entire  water 
column.  In  this  instance  the  circular  array  may  be  regarded  as  a  line 
array  with  a  circular  aperture  weighting  function.  The  result  of  this 
"weighting"  function  is  to  produce  an  azimuthal  beam  pattern  with  the 
beamwidth  of  a  linear  array  having  a  circular  weighting.  From  Refer¬ 
ence  C-2,  this  is 


circ 


je  x/L  =  1.13  B 


lin 


(C-10) 


Using  this  result  in  Equation  C-9  yields  the  error  estimate 

°-316(Bcir43/2 


(To 


Vi'Fy' 


(C-ll) 
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In  the  case  of  a  circular  array  with  a  baffle  to  eliminate  noise  inputs 
from  the  180°  sector  behind  the  array,  there  is  a  factor  of  2  increase  in 
the  effective  SNR  per  element  relative  to  omni-directional  noise.  As  this 
is  the  case  of  actual  interest  to  us,  we  write 


(C-12) 


For  convenience,  we  express  botho^  and  B  ..  in  terms  of  degrees 
rather  than  radians  by  dividing  both  by  180°  /n: 


(C-13) 


We  now  consider  a  specific  application.  Let  the  integration 
time  T  and  bandwidth  F  correspond  to  a  single  tone  with  TF  product  of  one. 
Than,  the  rms  error  becomes 


(C-14) 


-This  process  is  now  repeated  for  a  total  of  N  tones .  Assuming  that  these 

t 

tones  are  all  received  independently,  we  are  then  averaging  over  N  samples, 
and  the  rms  error  thus  decreases  as ^nT  Our  final  expression  for  the 
error  (in  degrees)  is  then  given  by 


(C-15) 


C-5 


Equation  C-14  is  plotted  in  Figure  C-l  with  cr  °Jtfy  as  a  function 

D 

of  B°circ  and  in  Figure  C-2  with  cr°  as  a  function  of  y,  parameterized  by  N, 
for  B°c£rc  equal  20°. 

Probably  the  most  significant  error  in  using  this  result  will  be 
due  to  the  omission  of  SNR  improvement  that  can  result  from  vertical  gain 
due  to  the  general  non-uniformity  of  the  noise  field  in  elevation.  This 
can  be  especially  significant  at  higher  frequencies.  It  should,  however, 
be  a  useful  limit,  to  the  accuracy  available  from  the  use  of  a  circular 
array  with  a  unidirectional  illumination  pattern. 
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ABSTRACT 

This  report  analyzes  the  sending  of  a  one-bit  message  over  a 
Rayleigh  fading  channel.  Transmission  of  N  bits,  selected  according  to 
a  frequency  hopped  pattern,  is  analyzed.  This  choice  is  based  on  con¬ 
siderations  of  the  channel  itself  and  a  presumed  covertness  requirement 
Performance  results,  in  the  form  of  detection  probability  curves,  are 
presented.  Some  new  (?)  results  for  calculating  and/or  inverting  the 
incomplete  Gamma  ratio  are  included. 
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SECTION  1 


INTRODUCTION 

The  purpose  of  this  report  is  to  present  an  approach  to  solving 
a  communications  problem,  explain  the  reasons  underlying  this  approach, 
describe  briefly  the  pertinent  mathematical  analyses  and  summarize  the 
results  as  a  set  of  representative  performance  curves.  We  consider  here 
the  problem  in  general;  application  to  a  specific  problem  follows  readily 
from  these  results. 

The  problem  consists  of  devising  and  analyzing  a  system  for 
transmitting  a  message.  This  message  is  equivalent  to  one  bit  of  infor¬ 
mation  -  e.g.,  it  might  be  an  "Alert",  "Ready",  or  "Acknowledge"  signal. 
The  message  is  seldom  sent.  When  it  is,  it  is  a  time  unknown  to  the 
receiver,  although  the  allowable  transmission  times  may  be  discrete 
instants,  e.g.,  every  minute.  In  order  to  receive  any  such  message,  the 
receiver  must  continually  "look"  for  the  alert  message.  We  then  have  a 
problem  analogous  to  the  classical  radar  detection  problem,  wherein 
performance  is  determined  by  the  probability  of  detection  versus  the 
probability  of  false-alarm. 

We  specify  two  significant  constraints  implicit  in  the  design 
of  the  transmission  system.  One  is  the  nature  of  the  communications 
channel.  For  this,  we  assume  a  fading  channel  with  considerable  multi- 
path  spread  and  frequency  smear.  Further,  the  nominal  time  delay  and 
frequency  shift  between  the  transmitter  and  receiver,  although  bounded, 
is  unknown.  Specifically,  this  might  be  an  underwater  acoustic  communica¬ 
tions  channel  with  moving  transmitter  and  receiver.  The  second  constraint 
concerns  the  possible  existence  of  unauthorized  receivers,  which  we  call 
interceptors.  This  means  that  we  wish  to  structure  the  message  such 


that  It  is  recognizable  to  the  receiver  but  not  (or,  more  exactly,  to  a 
much  lesser  extent)  to  the  interceptor. 

The  above  considerations  lead  to  a  redundant  message  structure. 
Although  the  message  has  information  content  of  only  one  bit,  we  picture 
it  as  consisting  of  a  collection  of  transmissions,  each  individual  trans¬ 
mission  being  a  pulsed  tone.  The  tones  are  scattered  throughout  a  region 
of  time  and  frequency  according  to  a  pattern  known  both  to  the  transmitter 
and  intended  receiver.  This  may  be  viewed  as  providing  time  and/or 
frequency  diversity  to  partially  compensate  for  the  fading  channel. 
Alternatively,  it  may  be  viewed  as  allowing  reliable  yet  covert  communica¬ 
tion  in  that  it  is  assumed  that  the  tone  structure  is  unknown  to  a  would-be 
interceptor.  Both  interpretations  are  correct  in  that  the  redundant  pulsed 
tones  may  be  transmitted  with  less  total  signal  energy  than  would  be 
required  for  a  single  pulsed  tone  for  the  same  probability  of  detection. 

First  we  discuss  the  tone  locations  in  time  and  frequency  and 
the  effect  that  various  patterns  have  on  interceptability.*  Then,  we 
present  briefly  the  type  of  processing  required  at  the  receiver.  We 
also  consider  here  the  multipath  spread  and  Doppler  smear  characteristics 
of  the  channel.  Finally,  we  examine  the  anticipated  system  performance, 
aa  determined  by  the  total  number  of  tones  comprising  the  message  and 
the  received  signal-to-noise  ratio  per  tone. 
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A  detailed  analysis  of  the  performance  realized  to  an  interceptor  shall 
be  presented  in  a  companion  report. 
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SECTION  2 


ALERT  MESSAGE  STRUCTURE 


2.1  THE  TIME /FREQUENCY  MATRIX 

Throughout  this  report  we  shall  utilize  a  two-dimensional 
(time  and  frequency)  representation  for  the  alert  message.  The  time 
dimension  will  be  divided  into  intervals  of  length  T,  referred  to  as 
"time-slices";  frequency  will  be  divided  into  intervals  of  size  Af, 
referred  to  as  "frequency-cells . "  For  simplicity,  we  shall  set  (Af)  T 
to  unity.  The  intersection  of  a  time  slice  and  a  frequency  cell  will  be 
referred  to  as  a  tone-location  or  a  tone.  The  intersection  of  N  continuous 
slices  and  M  cells  constitutes  a  "time-frequency  matrix"  of  size  N  x  M. 

An  example  of  such  a  matrix  is  shown  in  Figure  2-1.  We  will  denote  a  matrix 
of  this  form  by  S  and  a  tone  within  the  matrix  by  S(i,j).  This  definition 
of  S  is  thus  far  independent  of  absolute  time  and  frequency  values.  We 
now  restrict  our  attention  to  an  arbitrary  but  fixed  range  of  frequencies 
of  total  width  F.  Time  is  considered  to  be  an  infinite  continuum  quantized 


into  increments  of  size  T.  Thus,  each  time  quantum  defines  the  beginning 
of  a  different  time/frequency  matrix.  Where  necessary,  this  will  be 


denoted  as  where  t^  =  kT  represents  the  absolute  time  value  corresponding 


to  the  index  i  =  1. 


Note:  whenever  indices,  subscripts,  etc.  are  used,  the  first  term  will 
specify  time;  the  second,  frequency. 


M  FREQUENCY  CELLS 


TIME /FREQUENCY  MATRIX.  N=M=8 


2.2  BASIC  CONSIDERATIONS  . 

We  define  the  alert  message  in  terms  of  a  time/ frequency  matrix 
of  size  N  x  M,  where  N  and  M  are  to  be  chosen  as  part  of  the  actual  signal 
design.  The  message  itself  is  defined  as  a  pattern  of  tones  within  this 
matrix  S.  Each  time  slice  of  S  determines  a  transmission;  within  each 
slice,  the  particular  tones  (one  or  more)  dictated  by  the  pattern  are 
activated.  Thus,  the  requisite  transmitter  for  this  system  is  similar  to 
that  which  would  be  used  for  multi-tone  On/Off  keying  (MTOOK) .  For  purposes 
of  analysis,  the  time  origin  of  the  message  is  irrelevant.  The  ordering 
of  active  tones  throughout  the  time/frequency  matrix,  known  as  the  pattern 
or  key  to  the  message,  is  all  that  is  important. 

It  was  earlier  stated  that  the  channel  was  subject  to  fading. 

We  note  here  that  the  delay  and  frequency  spread  are  such  that,  in 
general,  coherent  combination  of  the  received  tones  is  impossible.  Thus, 
incoherent  combination  must  be  used.  Let  the  total  number  of  tones  in 
the  pattern  be  L.  As  long  as  the  key  is  known  to  both  the  transmitter 
and  receiver,  and  the  time/frequency  matrix  is  structured  so  that  each 
cell  is  subject  to  independent  Rayleigh  fading,  the  actual  location  of 
these  L  tones  within  the  N  x  M  matrix  does  not  matter.  That  is,  system 
performance  will  be  dictated  purely  by  the  total  number  of  tones  and  the 
total  transmitted  signal  energy.  The  assumed  independent  fading  is 
realistic,  especially  if  the  individual  tones  are  sufficiently  scattered 
in  time  and/or  frequency.  If  not,  if  there  is  correlated  fading  among 
the  tones  in  the  pattern,  full  diversity  effects  will  not  be  realized 
and  performance  will  be  somewhat  degraded  compared  to  the  idealized  case. 

Scattering  of  the  individual  message  tones  in  time  and 
frequency  is  in  basic  agreement  with  our  desire  to  make  the  signal 
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non -detectable  to  an  unauthorized  receiver.  Intuitively,  it  is  clear 
that  the  greater  the  time/ frequency  space  within  which  the  L  tones  are 
scattered  and  the  less  the  clustering  of  tones  within  that  space,  the 
less  detectable  the  message  will  be  to  someone  not  possessing  the  key. 

We  assume  that  the  pattern  chosen  is  sufficiently  diverse  to  ensure 
independent  fading  regardless  of  its  detailed  composition.  We  then  more 
finely  define  the  message  key  based  on  considerations  of  interceptability. 

We  do  not  intend  here  to  derive  expressions  for  intercept 
probabilities,  rather,  we  wish  merely  to  derive  a  general  conclusion 
concerning  the  form  of  the  pattern  such  that  the  probability  of  inter¬ 
cept  will  be  minimized.  Although  the  interceptor  does  not  know  the 
pattern  we  will  assume  that  it  has  perfect  knowledge  of  the  time/frequency 
matrix  -  i.e.,  it  knows  T,  Af  ,  N  and  M  and  is  synchronized  in  both  time 
and  frequency.*  Then,  for  each  time  slice  it  receives,  it  may  separately 
detect  the  energy  in  each  frequency  cell.  Alternatively  it  may  detect 
the  energy  contained  in  the  total  band. 

Whether  narrow  or  wide-band  detection  is  used ,  we  will  assume 
for  present  purposes  that  the  interceptor  examines  only  one  time  slice 
at  a  time,  each  slice  independent  of  all  others.  One  might  identify 
more  sophisticated  processing  schemes  for  the  interceptor  wherein  the 
entire  matrix  is  examined.  These  will  be  examined  in  a  companion  report. 
Further,  we  note  that  the  intercept  probability  for  many  schemes  of  this 
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form  will  be  directly  and  monotonically  related  to  the  intercept  probability 
for  one  time  slice.  Thus,  it  is  meaningful  to  structure  the  message  to 
minimize  the  basic  intercept  probability  P^.  for  one  slice. 

Our  initial  discussion  showed  the  detection  probability  to  be 
a  function  of  the  total  signal  energy  E  and  the  number  of  tones  L  in  the 
message.  This  fixes  the  energy  per  transmitted  tone  at  E/L.  Given  this 
constraint,  it  is  clear  that  the  least  interceptable  signal  is  one  con¬ 
taining  an  average  of  one  or  fewer  tones  per  time  slice.  To  see  this, 
consider  briefly  two  ways  in  which  an  unauthorized  receiver  might  process 
each  time  slice.  If  the  interceptor  employs  narrow-band  analysis,  any 
active  tone  will  be  equally  detectable  with  probability  P;  the  probability 
of  detecting  one  or  more  out  of  n  tones  will  be  of  the  form  1  -  (1  -  P)n 
which  increases  with  n.  For  wide-band  analysis,  increasing  the  number 
of  tones  increases  (linearly)  the  total  energy  in  the  band  -  once  again, 

P^.  increases  monotonically  with  the  number  of  tones.  Thus  our  conclusion.* 
The  general  nature  of  the  above  results  leads  to  the  somewhat 
heuristic  conclusion  that  the  alert  message  should  be  structured  such 
that  in  any  time  slice,  at  most  one  tone  is  active.  Just  as  examination 
of  one  time  slice  (all  the  tones  in  one  row)  of  the  matrix  was  postulated 
for  the  interceptor,  so  too  could  we  consider  examination  of  one  frequency 
cell  (all  the  tones  in  one  column).  This  leads  to  a  conclusion  which  is 
the  logical  dual  of  the  above  -  in  any  frequency  cell,  at  most  one  tone 
should  be  active.** 

*If  one  assumes  a  total  transmitted  energy  per  time  slice  regardless  of 
the  number  of  active  tones,  a  different  result  will  follow  (Reference  2). 
This  does  not  apply  to  our  problem. 

**Subject  to  practical  limitations  as  discussed  later. 
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Combining  these  criteria,  we  require  a  message  pattern  consisting 


of  at  most  one  active  tone  in  any  row  or  column  of  the  matrix.  This  is 
conceptually  equivalent  to  the  problem  of  placing  rooks  on  a  chessboard 
such  that  no  rook  can  capture  any  other.  Many  solutions  are  possible  - 
an  example  is  shown  in  Figure  2-2. 

The  above  result  is  certainly  not  surprising.  It  simply 
represents  the  greatest  extent  of  scattering  of  the  tones  in  the  allowed 
region.  It  leads  to  the  obvious  conclusion  that  either  dimension  of  the 
matrix  must  be  at  least  equal  to  the  number  of  tones  in  the  message, 

N  >  L  ;  M  >  L  (2-1) 

Having  an  inequality  in  Equation  (2-1)  is  equivalent  to  having  a  totally 
unutilized  time  slice  or  frequency  cell.  This  spreads  the  message  over 
an  even  greater  time/frequency  area.  Although  this  might  decrease 
intereeptauility  when  more  sophisticated  techniques  are  employed,  it 
also  increases  the  operational  requirements  for  the  intended  receiver. 

We  will  assume  that  our  message  is  sufficiently  diverse  and  that  we  do 
not  wish  to  make  the  message  length  or  bandwidth  any  larger  than  is 
necessary.  We  then  have 

L  =  N  =  M  (2-2) 

which  we  will  assume  to  be  true  in  the  rest  of  this  report.  The  symbol 
N  will  interchangeably  be  used  to  denote  the  number  of  tones  in  the 
message  or  the  time  or  frequency  dimension  of  S. 

It  should  be  noted  that  the  above  arguments  tacitly  assume 
that  a  message  structure  as  defined  by  Equation  (2-2)  or  (2-3)  will  not 
be  of  excessive  duration  or  bandwidth.  If  the  number  of  transmitted 
tones  required  to  yield  a  desired  level  of  performance  is  too  large,  then 
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more  than  one  active  tone  must  be  put  into  some  rows  or  columns  of  the 
matrix.  As  the  communications  system,  especially  for  some  of  the  channels 
under  consideration,  is  bandwidth  limited,  we  "double-up"  in  time.  That 
is,  we  hold  to  only  one  tone  per  time  slice  but  allow  more  than  one  tone 
per  frequency  cell.  Presumably,  the  message  might  be  spread  in  time 
across  n  matrices,  each  of  size  N  x  M.  Each  submatrix  would  be  as  dis¬ 
cussed  in  this  section,  the  composite  matrix  would  be  of  size  nN  x  M  and 
the  message  would  contain  nN  toties.  Our  formulations  for  system  performance 
will  not  be  affected  by  this;  one  might  wish  to  modify  the  expressions  for 
Pj  to  account  for  the  extra  detectability  introduced  by  this  modified 
structuring.  Under  most  conditions,  however,  enough  bandwidth  should  be 
available  to  accommodate  the  pattern  within  one  N  x  M  matrix,  and  such 
refinements  will  not  be  needed. 
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SECTION  3 


PROCESSING  THE  RECEIVED  SIGNAL 

In  this  section  we  describe  the  form  of  processing  used  at  the 
receiver  to  determine  whether  a  message  is  present  or  not.  The  implications 
of  the  time  and  frequency  spreading  and  shifting  associated  with  the 
channel  are  also  discussed. 

3.1  THE  RECEIVER 

The  receiver  constructs  time/ frequency  matrices  identical  to 
that  determining  the  transmitted  signal.  One  matrix  is  formed  every  T 
seconds.  The  individual  tone  cells  are  of  duration  T  and  width  Af  =  1/T 
as  discussed  earlier.  In  order  to  do  this,  input  samples  must  be 
accumulated  at  a  rate  M/T;  these  are  then  processed  via  an  FFT  to  produce 
M  frequency  cells  for  one  time  slice.  Each  time  (t^)  defines  a  new 
matrix  consisting  of  the  N  most  recent  time  slices. 

From  each  matrix,  the  tones  corresponding  to  the  predefined  key 
are  extracted  and  square-law  combined  to  obtain  a  decision  variable  which 
is  compared  to  a  threshold  level.  Th-is  represents  optimum  (maximum 
likelihood)  processing  for  independent  Rayleigh  fading  signals  in  the 
presence  of  white  Gaussian  noise.  The  processing  operations  are  illustrated 
in  Figure  3-1,  wherein  the  time/frequency  matrix  is  depicted  as  M  shift 
registers  in  parallel,  each  of  length  N.  The  post-detection  shift  register 
processing  produces  an  output  which  corresponds  to  the  cross -correlation 
of  the  message  key  with  the  detected  power  level  in  the  resolution  cells. 

Let  be  the  matrix  containing  all  N  received  message  tones. 
Although  Sj^  and  each  contain  N-l  of  the  N  transmitted  tones,  they 


3-1 


SELECT  I  _  I  THRESHOLD 


are  out  of  alignment  with  the  pattern.  The  particular  form  of  pattern 
specified  in  the  previous  section  is  such  that  any  displacement  from 
perfect  alignment  results  in  square-law  combining  the  envelopes  of  N 
resolution  cells,  none  of  which  contain  signal.  Thus,  the  tones  selected 
according  to  the  key  from  a  received  matrix  ideally  either  all  contain  a 
signal  (and  this  happens  for  only  one  time  t^) ,  or  are  all  purely  noise. 
Thus,  the  correlation  function  formed  by  the  receiver  is  ideally  a  delta 
function  when  a  message  is  present. 

In  actuality,  there  will  be  spillover  both  in  time  and  frequency 
which  will  broaden  the  correlation  function.  This  results  both  from  the 
spreading  and  the  shifting  of  the  individual  tones  comprising  the  message. 
Nonetheless,  the  receiver  responses  at  T  seconds  before  and  after  the 
peak  are  likely  to  be  small  compared  to  the  response  at  the  peak.  Thus, 
the  correlation  function,  although  not  an  impulse,  will  still  be  quite 
narrow.  The  significance  of  this  is  that,  in  order  to  have  a  reasonable 
probability  of  detection,  the  correlation  processing  depicted  in  Figure  3-1 
must  be  performed  every  time  slice. 

A  consequence  of  the  sharpness  of  the  receiver  correlation 
function  output  when  a  message  is  present  is  worthy  of  note.  Our  analysis 
will  assume  a  stationary  noise  environment.  Then,  when  no  signal  is 


present,  the  decision  variables  calculated  from  successive  matrices 
and  + ^  will  be  independent.  If  the  threshold  is  exceeded  at  time 


it  most  likely  will  not  be  exceeded  at  either  t^  ^  or  t^  +  ^.  That  is,  at 


V 


the  correlation  function  will  be  similar  in  shape  to  that  when  a 


signal  is  present  (it  will  have  a  sharp  peak)  and  will  be  indistinguishable 


from  that  due  to  a  true  message.  If,  however,  the  noise  environment  is 
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non-stationary ,  or  if  there  is  interference,  then  the  noise-only  outputs 
will  be  correlated  and  several  successive  time  instants  may  result  in 
threshold  exceedances.  Separately  considered,  these  would  each  be  false 
alarms.  Considered  as  a  whole,  they  form  too  broad  a  peak  for  the 
correlation  function  and  can  be  rejected  as  being  non-signal  in  origin. 
Thus,  by  a  small  amount  of  extra  processing,  one  can  improve  performance 
when  non-stationary  noise  or  interference  is  present.  In  effect,  one  is 
correlating  the  correlation  function  to  reject  a  certain  class  of  false 


alarm. 


3.2  EFFECTS  OF  THE  CHANNEL 

In  this  section  we  describe  qualitatively  the  changes  imparted 
to  the  individual  message  tones  by  the  channel.  The  receiver  processing, 
described  in  the  previous  section  and  analyzed  later,  assumes  these  to 
be  negligible.  Actually,  these  effects  can  be  significant,  and  the 
idealized  results  should  be  modified  to  account  for  these  disturbances. 

The  significant  impairments  to  the  message  due  to  the  channel  are  described 
in  Section  3.2.1  and  3.2.2.  In  Section  3.2.3  we  then  discuss  briefly 
how  to  modify  the  ideal  results  to  incorporate  these  effects  for  analysis 
of  an  actual  system. 

3.2.1  Spreading 

The  original  transmitted  tones  have  a  specified  duration  and 
frequency  width.  These  are  both  increased  due  to  passage  through  the 
channel.  This  results  in  spillover  of  signal  energy  into  adjacent  tone 
locations  in  the  received  time/ frequency  matrix.  As  the  matrix  is  sparsely 
inhabited,  this  is  not  as  significant  a  problem  as  it  might  be  in  a 
conventional  multitone  communications  system.  However,  as  the  received 
signal  energy  is  spread  over  a  greater  time /frequency  area,  less  signal 
energy  is  contained  within  the  nominal  tone  location  in  the  matrix. 

This  is  effectively  modeled  as  a  decrease  in  the  received  signal-to-noise 
level  per  tone. 

3.2.2  Shifting 

The  alert  message  is  defined  by  a  pattern.  This  pattern  is 
relative — it  is  made  absolute  by  choosing  an  origin  for  it  in  the  time/ 
frequency  plane.  The  message  is  received  with  this  origin  shifted  The 
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time  shift  is  caused  by  propagation  delay;  the  frequency  shift,  by  the 
effects  of  Doppler.  We  assume  that  the  amount  of  these  shifts  are  unknown, 
but  that  bounds  for  them  are  available.  The  receiver  must  then  be  able 
to  detect  the  message  within  some  rectangular  region  of  uncertainty  in 
time/ frequency  centered  about  a  nominal  origin. 

Since  the  signal  is  sent  at  a  time  unknown  to  the  receiver,  the 
nominal  origin  in  time  of  the  received  signal  is  unknown.  Thus,  shifting 
this  origin  in  time  has  no  effect  on  the  receiver  as  it  must  search 
through  all  time  for  the  message  anyway.  We  note  here  that  the  only  instance 
in  which  this  time  shifting  is  significant  is  when  the  propagation  path 
fluctuates  so  rapidly  that  different  time  slices  in  the  message  are 
delayed  by  different  amounts.  We  assume  that  this  does  not  happen---i.e. , 
there  is  not  "time  distortion."* 

In  order  to  avoid  searching  throughout  the  frequency,  as  well 
as  the  time  continuum,  it  is  assumed  that  the  frequency  origin  of  the 
message  is  fixed  and  known.  Then,  ideally,  the  receiver  need  look  only 
for  the  pattern  as  defined  by  this  nominal  frequency  value. 

In  actuality,  the  message  is  received  in  a  frequency  band 
related  to  the  nominal  band  by  a  shifting  up  or  down  in  frequency  and  a 
simultaneous  stretching  or  compressing.  If  the  relative  velocity  between 
the  transmitter  and  receiver  were  known,  then  the  Doppler  offset  (as  a 
function  of  frequency)  would  be  known  and  the  pattern  of  the  received 
message  tones  could  be  determined  from  the  nominal  pattern.  The 
receiver  would  then  merely  examine  the  time/ frequency  matrix  for  a 

* 

In  the  design  of  an  actual  system,  one  has  to  consider  the  change  in  path 
length  from  message  start  to  finish  caused  by  motion  of  the  transmitter 
relative  to  the  receiver.  This  may  not  be  so  readily  dismissed. 
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message  defined  by  this  distorted  pattern. 


We  assume  that  this  knowledge  is  not  available-— all  that  is 
known  is  a  maximum  possible  relative  velocity,  This  in  turn 

specifies  a  maximum  possible  frequency  shift.  Let  the  highest  transmitted 
frequency  be  .  Then,  the  maximum  shift  is  imparted  to  this  tone  and 
is  bounded  by 

SlAX  ^ 

where  c  is  the  propagation  velocity.  When  this  is  large  enough  to  cause 
the  received  tone  corresponding  to  to  lie  in  one  of  several  possible 

frequency  cells,  examination  in  the  receiver  of  multiple  Doppler  hypotheses 
is  necessary  in  order  to  ensure  detection  of  the  message. 

As  the  Doppler  shift  is  not  constant  across  the  band  of  frequencies 
comprising  the  message,  formation  of  these  multiple  hypotheses  is  not  a 
straightforward  procedure.  The  optimum  hypothesis  test  involves  the  use 
of  a  separate  sampling  frequency  in  the  receiver  processing  for  each 
postulated  value  of  (V/c).  The  processing  burden  and  memory  requirements 
(storage  of  NM  terms  for  each  hypothesis)  presumably  render  this  approach 
unfeasible.  A  cruder  but  simpler  method  is  to  compute  the  received  time 
and  frequency  of  each  tone  under  a  postulated  value  V/c  and  to  select 
the  nearest  time/ frequency  cell  from  the  nominal  matrix.  This  results 
in  sampling,  trans forming,  and  storing  the  incoming  signal  once  and  then 
.  performing  multiple  correlations  on  it.  Each  correlation  process  corresponds 
to  a  particular  hypothesized  distortion  of  the  nominal  message  pattern. 

The  number  of  such  correlations  required  may  be  equated  to  the 
number  of  cells  in  which  the  highest  tone  can  be  received.  From 
equation  (3-1),  this  is  roughly 

#  hypotheses  (2  AF/Af)  (3-2) 


More  exactly,  the  nominal  pattern  will  always  be  examined.  As  |Af  I 
increases,  an  equal  number  of  patterns  representing  positive  and  negative 
values  of  Af  will  be  required.  Then, 

Af  i 

#  hypotheses  =1+2  ^  +  —  (3-3) 

where  L*J  is  the  "greatest  integer  less  than  or  equal  to"  function. 

A  related  problem  due  to  shifting  is  lack  of  synchronization. 

As  the  message  appears  only  once,  there  is  no  chance  for  the  receiver  to 
adaptively  change  its  timing  and  sampling  rate  to  establish  sync.  Thus, 
a  received  tone  may  straddle  two  time  slices.  In  the  presence  of  Doppler 
it  may  also  straddle  two  frequency  cells.  This  would  be  true  even  if 
there  were  no  spreading  in  time  or  frequency.  Thus,  only  part  of  the 
available  signal  energy  is  used  when  a  particular  time /frequency  cell 
is  examined. 

This  could  be  partially  alleviated  by  finer  quantizing  in  time 

and/or  frequency - i.e.,  rather  than  rounding  a  postulated  received 

frequency  to  the  nearest  Af,  the  response  at  the  nearest  (Af/2),  etc., 
could  be  used.  The  additional  resolution  in  frequency  is  easily  obtained 
by  means  of  decimated  FFT  processing.  An  equivalent  increase  in  time 
resolution  can  be  provided  by  FFT  overlap.  Doubling  the  resolution  in 
each  dimension  increases  the  FFT  processing  by  slightly  more  than  a 
factor  of  four  and  the  required  buffer  space  by  a  factor  of  four. 

Whether  this  processing  burden  is  warranted  in  order  to  achieve  the 
corresponding  performance  gain  (less  than  2dB)  depends  on  the  performance/ 
cost  tradeoffs  for  the  particular  problem  at  hand.  Similarly,  even  finer 
resolution  could  be  obtained,  but  one  rapidly  reaches  the  point  of  vanishing 
returns. 
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3.2.3  Interpretation  of  Results 

In  Section  4  we  will  derive  results  applicable  to  a  non-dispersive 
fading  channel  and  a  receiver  which  is  synchronized  to  the  received 
message  tones.  To  interpret  these  results  in  a  more  realistic  setting, 
modifications,  as  indicated  in  the  discussion  of  spreading  and  shifting, 
have  to  be  made.  The  form  of  the  given  results  remains  valid,  but  the 
various  parameters  in  terms  of  which  these  results  are  formulated  need  to 
be  adjusted. 

One  necessary  modification  results  from  the  multiple  correlations 
required  to  test  the  various  possible  Doppler  hypotheses.  Let  the  number 
of  such  be  k.  We  stated  earlier  that  the  receiver  processed  the  incoming 
data  once  every  T  seconds.  This  is  still  true,  but  now,  k  rather  than 
1  is  the  number  of  decisions  made  every  T  seconds.  Thus,  the  effective 
time  interval  between  tests  is  T/k.  A  common  measure  of  performance  for 
a  detection  system  is  the  false-alarm  interval,  which  we  define  as  the 
time  interval  during  which  the  probability  of  having  one  or  more  false 
alarms  is  one  half.  This  is  a  function  of  the  false  alarm  probability 
per  decision  and  the  rate  at  which  the  decisions  are  made.  When  using 
the  performance  curves  given  in  Section  4,  it  should  be  remembered  that 

is  defined  per  decision,  not  per  time  slice.  The  false  alarm  interval 
is  then  determined  from  Table  4-1  by  using  the  value  of  T/k,  not  T. 

The  number  of  multiple  hypotheses  required  clearly  depends  on 
the  choice  of  tone  width.  However,  the  false  alarm  interval  is  essentially 

To  be  exact,  the  k  multiply  formed  decision  variables  are  not  independent. 
Thus,  use  of  T/k  is  a  worst  case  adjustment.  Dividing  T  by  only  one-half 
of  k  might  be  more  realistic. 
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independent  of  this.  To  see  this,  let  k  be  given  by  equation  (3-2).  As 
the  time/ frequency  matrix  is  constructed  with  Af  =  1/T,  we  alternatively 
have 

k  ~  (2AF)T 

and  the  effective  time  interval  between  decisions  is 


The  choice  of  basic  tone  size  thus  has  little  impact  on  the  system  operating 
characteristics  as  determined  by  the  false  alarm  interval.  It  might  seem 
that  the  number  of  multiple  hypotheses  would  be  a  measure  of  the  requisite 
complexity  of  the  receiver.  More  significant  is  the  rate  at  which  these 
decisions  must  be  made,  which  determines  the  required  processing  capability. 
From  equation  (3-4),  this  too  is  independent  of  the  basic  tone  structure. 

We  now  consider  spreading  and  lack  of  synchronization.  Although 
different  phenomena,  they  are  both  manifested  by  a  decrease  in  the  signal 
energy  in  the  time/f requency  cell  examined.  The  process  discussed  to 
counter  an  unknown  Doppler  shift  is  such  that  the  center  (in  time  and 
frequency)  of  a  received  message  tone  is  not  necessarily  centered  in  the 
selected  time/frequency  'ell.  Rather,  it  is  randomly  located  in  the 
cell,  uniformly  distributed  in  both  time  and  frequency. 

For  a  simple  analysis  of  the  magnitude  of  this  effect,  we 
consider  the  transmission  of  a  Gaussian  pulse.  We  assume  it  is  received 
without  distortion.  Using  matched  filtering  at  the  receiver,  the  output 
is  equivalent  to  the  sampling  of  a  bivariate  (both  time  and  frequency) 
Gaussian  function.  When  the  receiver  is  perfectly  synchronized ,  this 
function  is  sampled  at  its  peak.  In  our  case,  we  are  sampling  at  a 
point  randomly  located  in  a  time/frequency  cell  centered  about  the  peak 


of  the  Gaussian  function.  The  expected  value  of  this  sampled  point  can 
be  shown  to  be  down  (relative  to  the  peak)  by  2.05  dB.  By  halving  the 
uncertainty  in  both  time  and  frequency,  as  discussed  earlier,  this 
synchronization  loss  can  be  reduced  to  .61  dB. 

The  loss  in  signal  energy  due  to  spreading  cannot  be  so  readily 
estimated.  It  depends  on  both  the  form  and  degree  of  the  spreading. 

Given  a  specific  model  for  the  spreading,  the  shape  of  the  received 
waveform  could  be  determined.  This  is  done  by  convolving  the  pulse 
(time  and  frequency  representations  thereof)  with  the  channel  spreading 
functions.  As  the  signal  is  stretched  out  (both  in  time  and  frequency), 
the  peak  value  must  decrease  in  order  to  maintain  the  original  level  of 
signal  energy.  This  represents  the  spreading  loss  in  the  received  signal. 

Combining  the  two  effects  yields  an  estimate  of  the  total  required 
adjustment  to  the  received  SNR  per  tone  to  account  for  shifting  and 
spreading  losses.  These  tend  to  mollify  each  other---as  the  received 
pulse  is  spread  out,  the  loss  incurred  by  sampling  off  the  peak  decreases. 
Thus,  the  net  effect  will  normally  be  less  than  that  which  would  be 
obtained  by  considering  the  two  effects  separately  and  independently 


combining  them. 


SECTION  4 


PERFORMANCE 

In  this  section  we  discuss  the  statistics  of  the  decision 
variable  formed  by  the  receiver  under  the  hypotheses  of  no  signal  present 
and  signal  present.  Two  types  of  threshold,  fixed  and  adaptive,  against 
which  the  decision  variable  is  to  be  compared,  are  discussed.  Probabilities 
of  false  alarm  and  detection  are  defined  for  both  cases.  The  mathematical 
formulations  are  briefly  discussed  and  results  are  presented  for  the  fixed 
threshold  case. 

As  discussed  earlier,  an  idealized  situation  is  assumed.  The 
channel,  although  fading,  is  non-dispersive ;  the  effects  of  Doppler  are 
loown  and  removed  at  the  receiver  and  the  receiver  is  synchronized  to  the 
message.  To  relate  these  results  to  a  particular  actual  example,  we  modify 
them  as  described  in  Section  3.2.3. 


4.1  INTRODUCTION 


Out  of  the  received  time/frequency  matrix,  N  complex  variables 

are  selected  according  to  the  pre-defined  tone  pattern.  Denote  these 

as  {u^,  1-i-N}.  We  consider  two  idealized  hypotheses:  1)  The  {u^}  are 

all  purely  noise;  2)  The  {u^}  are  all  signal  plus  noise. 

The  noise  terms  are  assumed  to  be  independent,  zero-mean,  complex 

2 

normal  random  variables  of  variance  o  .  We  are  considering  only  fading 

channels,  hence,  we  assume  the  signal  terms  to  also  be  independent,  zero- 

2 

mean,  complex  normal  random  variables  of  variance  og  .  As  the  individual 
tones  are  separated  in  time  and  frequency,  the  assumptions  of  independence 
are  realistic.  We  then  write  the  two  hypotheses  as: 

!•  ui  =  xi  +  jyi  •  xi»  y±  e  N(°.°2) 

2-  ui  =  xi  +  jyi  ;  xt,  yt  e  n(o,o2  +  os2) 

2 

where  the  notation  "x_^eN(p,a  )"  is  read  as  "x^  is  a  normal  random  variable 

o 

of  mean  p  and  variance  o  . " 

Detecting  a  signal  under  these  conditions  is  a  classic  problem. 
The  optimum  receiver  may  be  shown  (Ref.  7)  to  be  one  in  which  a  decision 
variable,  formed  by  summing  the  squared  magnitudes  of  the  tone  responses, 
is  compared  against  a  threshold  determined  by  the  variance  of  the  noise. 

We  denote  this  decision  variable  by  U, 

A  9  o 

U  =2  x.  +  y.2 
i=l 


As  the  {x^}  and  the  (y^)  are  identically  distributed  and  mutually 
independent,  we  may  equivalently  write 


U 


(4-1) 
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Note  that  U  is  essentially*  a  chi-squared  random  variable  with  2N  degrees 
of  freedom. 


The  signal-to-noise  ratio  is  defined  by 

2 


Y  = 


(4-2) 


and 


a  = 


1+Y 


2,  2 

a  +a 

s 


(4-3) 


The  two  hypotheses  defining  the  distribution  of  the  {x^}  in 
equation  (4-1)  may  then  be  written  as 


e  N(0,o  ) 


,  noise  only 


x^  e  N(0,o^/a)  ,  signal 


(4-4) 


The  test  to  determine  whether  or  not  an  alert  message  is  present 
consists  of  comparing  the  decision  variable  U  against  a  threshold,  which 
we  denote  as  W.  When  U  exceeds  the  threshold  but  a  signal  is  not  present, 
we  have  a  false  alarm;  inversely,  when  U  does  not  exceed  the  threshold 
but  a  signal  is  present,  we  have  a  missed  detection.  The  probabilities 
for  each  of  these  forms  of  error  are: 

P_.  =  P  {U>W  I  no  signal}  (4-5) 

FA  r  » 

PMD  =  Pr  *U<W  |  si8na1^  (4-6) 

The  probability  of  a  successful  detection  is 

pD  =  1  -  PMD  =  pr  (U>W|  signal}  (4-7) 

The  particular  probability  distributions  to  be  used  here  depend 
upon  the  form  of  thresholding  used,  as  described  in  the  next  two  sections. 

*Strictly  speaking,  a  chi-squared  random  variable  is  formed  by  summing 
the  squares  of  unit  variance  normal  random  variables. 


The  traditional  design  approach  in  a  detection  problem  is  to 
first  determine  the  maximum  false-alarm  rate  at  which  one  is  willing  to 


operate.  Given  this  value  of  P  ,  equation  (4-5)  is  inverted  to  deter¬ 
mine  the  required  threshold  W.  Using  this  threshold  in  equation  (4-7) , 
one  then  calculates  the  detection  probability  as  a  function  of  the 
received  signal-to-noise  ratio. 

For  this  particular  problem,  the  number  of  tones,  N,  is  not 
specified.  Thus,  our  objective  is  to  create,  for  a  given  P  ,  a  set  of 
curves  parametrized  by  N  depicting  detection  probability  as  a  function 
of  signal-to-noise.  Given  a  desired  level  of  performance,  described  by 
the  probability  of  detection  realized  at  a  certain  signal-to-noise  ratio 
and  false-alarm  rate,  the  appropriate  number  of  message  tones  may  then 
be  determined.  In  actual  practice,  one  would  also  have  to  determine 
the  intercept  probabilities  corresponding  to  the  chosen  parameters  in 
order  to  fully  evaluate  a  particular  signal  design. 
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4.2  THE  USE  OF  A  FIXED  THRESHOLD 


4.2.1  Formulation 

In  this  section  we  examine  the  use  of  a  fixed  value  for  the 

2 

threshold  W.  This  threshold  is  set  to  a  multiple  of  a  ,  the  variance 

of  the  noise,  which  implies  having  precise  knowledge  of  this  parameter. 

This  is  an  idealization,  yielding  optimum  performance.  In  the  next 

section  we  will  consider  an  adaptive  form  of  thresholding.  This 

2 

effectively  derives  an  estimate  of  a  from  measurements  of  background 
noise  and  represents  a  more  realistic  thresholding  technique. 

First  we  consider  in  greater  detail  the  decision  variable  U 
under  the  hypothesis  that  no  signal  is  present.  From  equation  (4-1) 
we  have 

2N  - 

U  =  x^  ,  xi  e  N(0,a  ) 
i=»l 


This  may  be  converted  to 


„  2N 

t  =  ~2  <X,-A0  ,  {x Jo)  e  N(0,1) 


1=1 

which  is  a  chi-squared  random  variable  with  2N  degrees  of  freedom.  The 
density  function  is 

N-l  -t/2 

ft  (t)  - 


2"  r(N) 


where  T(N)  is  the  gamma  function 


dt 


f  N-l  -t 
r(N)  =  J  t  e 

o 

r(N)  ■  (N-l)!  ,  N  an  integer 


(4-8) 
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Alternatively ,  the  decision  variable  U  may  be  expressed  in  terms  of  the 
random  variable  v, 


V  "  2 


U 

2  a2 


(4-9) 


v  is  said  to  have  a  gamma  distribution,  its  density  is  given  by 

(4-10) 


f  (v) 
v 


VN-1  e-v 


T(N) 


We  shall  use  this  form  of  the  decision  statistic. 

Next,  we  consider  the  decision  variable  under  the  alternative 
hypothesis  (signal  present) .  We  now  have 


U  =5^  x.2  ,  x.  e  N(0,cj2/a) 
i=l 


o 

As  before,  we  form  v  =  U/2a  •  The  presence  of  a  signal 
necessitates  an  additional  scale  factor  to  convert  the  decision  variable 
to  one  expressed  in  terms  of  unit  variance  normals."  Hence, 

v  =  £  [lE  (^  Va)2]  =  a  V'  (4_11) 


The  density  of  v  is  trivially  related  to  that  of  vr,  further,  the 
density  of  v'  is  the  same  as  that  for  -v  under  the  noise  only  hypothesis 
as  given  by  equation  (4-10) 

Thus, 


fv(v) 


signal, a 


af  (av) 
v 


(4-12) 


noise  only 


defines  the  probability  density  function  of  the  statistic  v  when  there 
is  a  signal  present. 


The  ideal  fixed  threshold  W  is  a  multiple  of  the  noise  variance, 

W  =  2Ao2  (4-13) 

Comparison  of  U  against  W  is  then  equivalent  to  the  test 

v  l  A  (4-14) 

with  v  defined  by  equations  (4-9)  and  (4-11) .  We  will  use  this  formula¬ 
tion  and  refer  for  simplicity  to  A  itself  as  the  threshold. 

Inserting  into  equations  (4-5)  and  (4-7),  we  get 


FA 


=  J  fy(v)  dv 


noise 


PD-  /’W'l 


dv 


signal 


By  (4-12),  the  latter  may  also  be  expressed  in  terms  of  the  density  for 
the  noise-only  hypothesis.  Dropping  the  designator,  the  expressions 
we  need  to  evaluate  become 


FA 


=  /  fv(v)  dv 


(4-15) 


f  Vv) 


dv 


(4-16) 


oA 


with  fv(v)  defined  in  equation  (4-10) . 

Although  the  formulations  presented  for  P  ^  and  are  con¬ 
ceptually  identical,  our  uses  of  them  are  not.  Equation  (4-15)  must  be 
inverted,  i.e. ,  given  Pp,^,  one  wishes  to  solve  for  A.  Then,  this  value 
Is  used  directly  in  equation  (4-16)  to  evaluate  P^. 

We  first  discuss  direct  evaluation  of  the  integrals.  This  is 
described  in  terms  of  P^,  although  it  also  applies  to  Pp.  Then,  inver¬ 
sion  of  the  integral  is  considered. 
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4.2.2  Solution 


Two  approaches  are  available  for  integration  of  the  density 
given  in  equation  (4-10) .  It  may  be  evaluated  as  a  series  consisting  of 
N  terms.  Specifically,  expressing  r(N)  as  (N-l)!,  it  is  trivial  to  show 


that 


N-l 

■feyr 


dv  =  -i 


-v 


i-0 


(4-17) 


Although  the  above  is  an  exact  form,  it  is  often  too  cumber¬ 
some  for  evaluation.  Examination  of  it,  however,  yields  some  useful 
insight.  The  summation  consists  of  the  first  N  terms  of  the  series 

representation  of  ev;  it  is  often  denoted  as  e  (v) .  In  this  form, 

N-l 

we  have 

PFA  '  «'#  Vl  <»>  (4-18> 

Examining  equation  (4-18),  we  see  that  for  a  fixed  A,  as  N 

goes  to  infinity,  P  approaches  unity.  In  retrospect,  this  is  not 
r  A 

surprising,  as,  in  order  to  maintain  a  given  Pj.A»  we  expect  A  to  grow 
with  N.  Actually,  it  is  mathematically  more  meaningful  (see  Appendix  A) 
to  express  A  in  terms  of  (N-l).  We  will  define  a  "normalized"  threshold 
8  by 

8  -  (4-19) 


Then,  we  may  write 

eN_1(8(N-l)) 

PFA  exp  (8 (N-l)) 

From  equation  (6.5.34) 

of  Ref.  4, 

0  , 

8>1 

lim 

(N-l) 

8  FA  '  112  ’ 

1  . 

8=1 

0s8<l 
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i 

4 


i 


Calculation  of  values  of  A  shows  that  8>1.  Hence,  as  the  message  length 
grows  infinitely  large,  the  false-alarm  probability  becomes  vanishingly 
small.  Thus,  we  can  expect  our  system  performance  to  improve  as  N 
increases . 


For  small  values  of  N,  evaluation  of  equation  (4-18)  is  feasible. 

As  N  (and  consequently  A)  grows  large,  the  magnitudes  of  the  terms 

become  increasingly  difficult  to  handle.  As  an  example,  values  of  N  -  32 

and  A  =  66.39  yield  (Ref.  5)  a  false-alarm  probability  of  10  Use  of 

-29 

equation  (4-18)  would  require  first  forming  the  terms  exp  (-A)  ~  10  ; 

I  O  O 

e__  ,  (A)  =  10  in  order  to  conclude  that  P_.  was  10  .  N  =  32  is 

N-i.  r  A 

hardly  to  be  considered  large,  in  fact,  we  are  interested  in  values  such 
as  N  =  100.  Then,  the  exponents  involved  in  (4-18)  are  such  that  main¬ 
taining  numerical  accuracy  can  become  a  serious  problem.  We  further 
note  that  inversions  of  (4-18)  is  clearly  extremely  awkward. 

Returning  to  equation  (4-15),  we  present  an  alternative  solution. 
Equation  (4-8)  defines  the  (complete)  gamma  function.  Similarly,  one 
has  the  incomplete  gamma  function  and  its  complement; 


y(N,x)  4  J*  tN  1  e  t  dt 
o 

T(N,x)  4  J  tN  1  e  dt 
x 

Y(N,x)  +  r(N,x)  =  TOO 


Then,  one  may  write 


p  =  r (n, A) 
FA  T(N) 


which  is  known  as  an  incomplete  gamma  ratio. 


(4-21) 


(4-22) 
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The  gamma  functions,  complete  and  incomplete,  have  been  rather 
thoroughly  studied,  and  many  approximations  or  expansions  are  tabulated 


(e.g.,  Refs.  4  and  6).  Surprisingly,  much  less  work  appears  to  have 
been  done  on  the  incomplete  gamma  ratio.  In  Appendix  A,  we  present 
what  we  consider  to  be  some  new  results  on  this  problem.  An  approximate 
solution  is  derived,  valid  for  calculation  of  the  false-alarm  probability, 
which  is  readily  invertable  for  determination  of  the  threshold  A.  The 
accuracy  increases  as  N  grows  large  or  as  the  desired  P  grows  small,  which 
are  the  two  conditions  that  aggravate  use  of  the  direct  summation  approach. 

This  approach  unfortunately  does  not  hold  for  calculation  of  the 
detection  probability,  except  for  very  small  signal-to-noise  ratios,  as 
the  asymptotic  expansion  upon  which  it  is  based  is  invalid  for  values 
of  the  argument  (aA)  that  appear  when  one  writes  equation  (4-16)  in 
the  form  of  an  incomplete  gamma  ratio.  Thus,  although  the  more  difficult 
problem  of  inversion  of  the  false-alarm  integral  is  solved  rather  easily, 
calculation  of  the  detection  probability  integral  must  be  done  in  a 
manner  similar  to  that  given  in  (4-18).  Actually,  a  modified  form  of 
this,  given  in  Ref.  5,  is  used.  This  is  also  explained  in  Appendix  A. 

4.2.3  Results 

In  the  next  few  pages,  we  present  families  of  performance 
curves  (Figures  4-1  through  4-7)  for  a  system  using  a  fixed  thresholding 
technique.  The  curves  are  presented  first  by  the  operating  level  at 
which  the  false-alarm  probability  is  set  and  then  by  different  values 
of  N,  the  message  length.  is  the  detection  probability  in  %  and  E/nc 
is  the  signal-to-noise  ratio  per  tone  in  dB. 
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FIGURE  4-5 


p 


A  simple  interpretation  of  is  given  in  Table  4-1.  Here, 
the  false-alarm  interval  (denoted  Tp^)  is  defined  as  the  length  of  time 
over  which  the  probability  of  at  least  one  false-alarm  occurring  is  one- 
half.  The  values  of  Tp^  given  are  normalized  to  a  tone  of  length  T  =  1  ms. 
(Note  that  although  the  overall  message  length  is  NT,  the  receiver  described 
in  Section  3.1  makes  a  decision  every  T  seconds.)  Thus,  for  instance,  if 
the  tone  length  is  20  ms  and  the  operating  level  of  the  system  is  set  for 
PpA  =  10  a  false-alarm  will  occur,  on  the  average,  every  3.8  hours. 
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TABLE  4-1 

FALSE  ALARM  INTERVALS 


t  =  False  Alarm  Interval,  normalized  for  T  =  1  ms 


p 

tfa 

3 

.69  seconds 

4 

6.9 

5 

1.15  minutes 

6 

11.5 

7 

1.91  hours 

8 

19.1 

9 

7.99  days 
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4.3  THE  USE  OF  A  VARIABLE  THRESHOLD 


4.3.1  Formulation 

In  this  section,  we  examine  the  use  of  a  variable  or  adaptive 
thresholding  technique.  Just  as  the  decision  variable  U  is  formed  by 
summing  the  squared  magnitudes  of  N  tones  which  may  or  may  not  contain  a 
signal,  so  too  is  the  threshold  W  formed  by  summing  the  squared  magnitudes 
of  N  tones  believed  to  be  purely  noise.  These  tones  may  be  derived  from 
a  frequency  band  wherein  no  signal  is  ever  expected,  alternatively,  they 
may  be  from  past  decisions  which  were  determined  to  represent  the  noise- 
only  hypothesis.  Actually,  one  would  create  the  noise  reference  from 
more  than  only  N  samples  in  order  to  obtain  a  better  estimate  of  the 
variance  of  the  noise  process.  Such  a  more  general  form  is  presented, 
the  implications  of  which  are  discussed  later. 

The  statistic  U  is  as  before, 

2 

x^  e  N(0,o  )  ,  noise 

2 

x.^  t  N(0,a  /a),  signal 


2N 


u=E 


i=l 


x. 

l 


The  threshold  W  is  given  by  the  average  over  v  intervals  of  length  N, 

V2N 


H  1  2 
w  =  v  £  wi 


(4-23) 


i=l 


where,  for  simplicity,  v  is  assumed  integral.  The  test  consists  of 
comparing  U  to  some  multiple  of  W,  specifically 


or 


U  J 


SW 


F  4  £  l  g 
W 


(4-24) 


Using  this  form,  we  will  refer  to  the  scalar  3  as  the  threshold  and  F 


as  the  decision  variable. 


The  decision  variable  may  be  expressed  as 


h  (u/°2> 

1 

2N 

2N 

Z 

(x±lo)2 

2 _  /  TJ  /  /-r  2  \ 

i=l 

2N  v*'0  ' 

”1“ 

v2N 

v2N 

I 

(w./a)2 

i 

Under  the  noise  only  hypothesis,  this  is  the  ratio  of  two  chi-squared 

variates,  each  normalized  by  its  number  of  degrees  of  freedom.  Such  a 

random  variable  is  known  as  an  F-distribution*  with  2N  and  v2N  degrees 

of  freedom.  Its  density  function  is  given  by 

f  (V)  =  (2N)1*  (u2N)vN  .  _ - 

FW  B(N,vN)  (v2N  +  2Nv) 

where  B  (m,n)  is  the  (complete)  beta  function  defined  by 

B(m,n)  ■  |  tm  ^  (l-t)n  *  dt  (4-25) 


We  thus  have 


PFA=  VV>dv 


(4-25) 


(4-26) 


We  derive  the  detection  probability  in  a  similar  manner.  The 
presence  of  the  signal  necessitates  a  different  normalization,  viz. 

r  k  1  (Ui,/‘,2)  ,ir. 
ii)  ("/o2) 

The  density  for  F'  is  the  same  as  that  for  F  under  the  noise  only 
hypothesis.  As  was  the  case  for  a  fixed  threshold,  has  a  formulation 
identical  to  P^, 


PD  e  J  fF  <v>  dv 


(4-27) 


More  specifically,  a  doubly  central  F  to  denote  that  each  of  the  chi- 
squareds  are  central  (formed  from  normal  r.v.s.  of  zero-mean). 


4.3.2  Solution 


It  can  be  shown  (Ref.  3c)  that  an  F-distribution  can  be 
expressed  as  a  beta  distribution.  Let  F  be  a  random  variable  having  an 
F-distribution  with  2N  and  v2N  degrees  of  freedom;  let  t  be  a  random 


variable  having  a  beta  distribution  with  parameters  N  and  vN.  Then, 


Pr  [F  <  K]  =  Pr  [t  <  K'] 


(4-28) 


where 


K'  =  ^  (4-29) 

The  probability  that  t<K^ is  given  by  an  incomplete  beta  ratio,  i.e. , 

Pr  [t  <  K']  =  IR.  (N,vN)  \ 


where 


B  (m,n) 

V"*10  =  TOa.nT 


and  is  the  incomplete  beta  function 
Bx  (m,n)  =  (*  tm'1(l-t)n_1d 


The  incomplete  beta  ratio  has  a  useful  symmetry  property 

Ix  (m,n)  =  1  -  Ij_x  (n,m) 

Using  these  identities,  we  may  then  write 
PFA  =  pf  tF  i  B]  =  1  -  Pr  [F  <  8]  \ 

PFA  "  1  "  V  <N’vN>  ( 


-  Ip„  (vN,N) 


with  6’  =  -£■ 


p  fill  _  V 

v+B  ’  v+B 


1-3’ 


and  similarly 
PD  =  Pr  [f  >  og] 

■  i- 1  «.,)•  <"•'«  ■  ‘w 


(vN,N) 


»lth 


1  -  (aB)’ 


(4-3C) 


|X  tm'1(l-t)n-1dt  •  (4-31) 


(4-32) 


(4-33) 


(4-34) 


,nd  m  u-B'f+o  T'  •  (ae)"  "  e,'-h.a-e"') 


The  formulations  involved  for  analyzing  performance  using  a 
variable  threshold  are  rather  similar  to  those  for  a  fixed  threshold, 
except  that  they  involve  beta  rather  than  gamma  functions.  Unfortunately, 
much  less  investigation  into  the  properties  of  the  beta  functions  appears 

to  have  been  done - at  least,  no  relations  resulting  in  a  feasible 

inversion  of  (4-33)  have  been  found. 

One  could  use  expression  (4-31)  for  the  incomplete  beta  function, 
expand  the  integrand  into  n  terms  and  integrate  each  directly.  This  is 
equivalent  to  summing  n  binomial  probabilities.  From  this,  the  incomplete 
beta  ratio  could  be  evaluated.  As  shown  in  the  next  section,  the 
computational  effort  involved  here  does  not  appear  to  be  justified. 

4.3.3  Discussion 

It  is  well  known  that,  as  the  variable  threshold  is  "noisier" 
than  an  ideal  fixed  one,  such  a  system  would  have  to  operate  at  a  higher 
threshold  setting  in  order  to  operate  at  the  same  false-alarm  rate  as  a 
fixed  system.  Thus,  one  would  expect  a  lower  detection  probability  and 
consequently,  poorer  performance  from  a  variable  system.  (This,  of 
course,  supposes  that  the  fixed  thresholding  technique  always  has  perfect 
knowledge  of  the  variance  of  the  noise  process.)  Further,  we  suspect 
that  forming  the  threshold  by  averaging  over  a  greater  number  of  observa¬ 
tions  will  provide  a  smoother  estimate  of  the  noise  and  hence  a  performance 
approaching  that  of  the  fixed  system.  It  is  easy  to  verify  these  intui¬ 
tions. 

We  earlier  discussed  the  fact  that  forming  the  threshold  over 
vN  observations  and  normalizing  by  1/v  allowed  a  greater  history  of  the 
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noise  process  to  be  used.  We  now  examine  the  behavior  of  such  a  threshold 

as  v  grows  large.  As  defined  in  equation  (4-23),  W  may  be  written  as 

u  c2  2 
“  '  —  *V2S 

i.e.,  a  chi-squared  with  (v2N)  degrees  of  freedome  multiplies  by  the 


scaling  factor  (o  /v) .  Since,  in  general, 

E  [/J  -  . 

var  h2,]  "  2“ 


we  have 


E  [w]  =  2No‘i 


[w]  = 


Letting  v  approach  infinity. 


^  E[tf]  *  2N°2 


[W]  =  0 


(4-35) 


Thus,  as  the  sample  size  is  increased,  the  "variable"  threshold  tends 
towards  a  fixed  threshold  as  discussed  in  Section  4.2.  Specifically, 

note  that  the  test  defined  in  equation  (4-24)  then  becomes 

?  2 
U  >  2Nga 

which,  other  than  a  minor  difference  in  notation,  is  identical  to  that 
presented  for  the  fixed  thresholding  technique. 

We  can  also  examine  the  system  performance  as  v-*».  For 
simplicity,  we  fix  N  =  1.  Then,  I  (v,l)  =  xV  and  the  equation  for  P  ., 

X  x*  A 

in  this  simple  case,  is  readily  inverted  for  g".  The  quantity  (ag)" 
may  be  expressed  in  terms  of  this  via  (4-34) ,  thereby  allowing  deter¬ 


mination  of  Pp.  For  large  v,g"  is  close  to  unity,  specifically. 


g"  -  1  -  e 


-In  Pv 


. . . 


mmmmm 


As  v«ot  it  is  then  easy  to  show  that  ^ 

PD  '  <PFA>“  '  <PFA>1+T  (4-36> 

which  is  precisely  the  result  one  would  get  for  a  fixed  threshold  analysis. 
This  same  approach  can  be  used  to  show  inferior  performance  for  finite 
values  of  v 

We  therefore  conclude  that  although  the  variable  thresholding 
method  will  produce  performance  inferior  to  a  fixed  threshold  for  small 
v,  in  actual  practice,  this  situation  will  not  occur.  Neither  will  an 
ideal  fixed  threshold  occur.  Rather,  an  implementable  system  will  have 
an  adaptive  thresholding  scheme  with  a  sufficiently  large  choice  of  v  to 
render  it  equivalent  to  a  fixed  system,  with  its  performance  then  as 
given  in  Section  4.2.3. 
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APPENDIX  A 


APPROXIMATIONS  TO  THE  INCOMPLETE  GAMMA  RATIO 

The  incomplete  gamma  ratio  is  defined  as  the  ratio  of  the 
complementary  incomplete  gamma  function  to  the  complete.  No  nomenclature 
is  established  for  this  quantity;  we  shall  denote  the  function  by  P(N,A), 
i.e. , 

P(N ,  A)  =  (A 

and  a  particular  evaluation  thereof  by  simply  P. 

We  consider  here  the  solution  of  two  related  problems.  The 
first  is  the  direct-given  N  and  A,  we  wish  to  find  P.  The  second  is  an 

inversion - given  P  and  N,  we  wish  to  find  A.  Approximations  for  these 

two  problems,  for  certain  ranges  of  values  of  N  and  A,  will  be  obtained 
appropriate  to  the  discussions  in  Section  4.2.2. 

We  first  define  a  normalized  threshold. 


The  choice  of  (N-l)  rather  N  as  a  normalizing  factor  is  not  logically 
obvious,  but  is  dictated  rather  by  the  formulations  which  follow.  We 
may  then  write  the  ratio  as  P(N,S). 

The  first  approximation  we  derive  is  valid  only  for  |3>1.  This 
is  of  interest  in  the  calculation  of  false-alarm  probabilities  and  has 
the  important  advantage  of  being  readily  inverted.  We  approximate 
P(N,A)  by  forming  the  quotient  of  asymptotic  expansions  for  each  of  the 
functions  appearing  in  (A-l) .  Although  these  expansions  are  separately 
known,  we  have  not  seen  this  combination  nor  the  useful  result  it  produces 


From  equation  9.5.4  of  Ref.  6,  we  have  the  asymptotic  expansion 


T(N , A) 


-A  .N 
e  A 

[A-(N-l)] 


[■ 


+  + » [-r^\]  <*-» 

LA-(N-l)]^  [a-(N-1)]J  L  [a-(N-1)]4-J_ 


This  is  valid  as  N  and  A  each  grow  large  provided  that  the  quantity 


z 


a  M 

'  A- (N-l) 


we  have 


goes  to  zero  from  the  positive  side. 


From  equation  (A-2) , 


[A-  (N-l)]  =  (N-l) (3-1) 


(A-4) 


Then,  z  =  y^j===  goes  to  zero  provided  that  3  does  not  go  to  unity. 

In  actuality,  3  does  approach  unity  and  for  values  of  $—1  the  approximation 
we  derive  will  be  incorrect.  However,  for  values  of  P  and  N  of  interest  to 
us,  (A-3)  is  seen  to  be  sufficiently  accurate.  Note  that  if  A< (N-l)  ,  as  would 
be  the  case  in  evaluation  of  detection  probabilities,  the  parameter  z  is  nega¬ 
tive  and  the  asymptotic  form  is  invalid.  Equation  (A-3)  may  be  re-written  as 


T(N,A) 


-A  ,N 
e  A 

(N-l) (3-1) 


1  - 


(N-l) (3-1) ' 


(N-l) 2 (3-1) 3 


+ 


0  - 

_(N- 


i)2(3-D4JJ 


(A-5) 


For  brevity,  we  will  express  the  term  in  brackets  as  [l  -  etc.]  .  We  will 
assume  some  if  not  all  of  the  "etc."  terms  negligible  in  the  following 
analysis.  Clearly,  if  the  terms  we  ignore  become  the  dominant 

ones . 


The  expansion  we  use  for  the  denominator  of  (A-l)  is  the  familiar 
Stirling  approximation: 

T(N)  =  /2?  (N-l)(N-1/2)  e"(N-1)  (1  +  small)  (A-6) 

The  terms  designated  "small"  are  tabulated,  we  will  maintain  only  the 
leading  term  (even  this  could  be  rather  safely  ignored)  which  is  1/12 (N-l). 


Combining  equations  (A-5)  and  (A-6) ,  we  now  have 


D  (N-l)(N_1/2)  /ST  (N-l) (3-1)  ~  (1  +  small) 


-(N-l) (3-1)  N  *r 

P(N,3)  s  - - E - (<v,X) 

(N-l) (3-1) 


where  "('''l)1’  denotes  the  ratio  of  the  last  two  terms.  This  is  further 
simplified  by  the  substitution 


e“  -  B  e01-1*  lnS 


Then, 


P(N,3) 


-  1  _3 _  -(N-l)  [(3-1) -In  3] 

~~  /nr  (N-i)  (3-d 


When  N  and  3  are  known,  this  may  be  readily  solved  for  P.  As  many  terms 
as  are  deemed  appropriate  may  be  maintained  in  (vl) . 


To  invert  equation  (A-7) ,  we  take  natural  logarithms  and  solve 


for  3. 


In  P  -  -(N-l)  [(3-1) -ln3]  +  ln3  -In (3-1) 

+  1/2  In  N  -  In  (N-l)  -1/2  In  (2tt)  +  In  (^1) 

The  term  (^1)  depends  upon  3;  to  get  an  initial  estimate,  we  ignore  it. 

To  simplify,  we  express  the  above  equation  in  terms  of  only  3  or  (3-1). 

That  is,  we  utilize  one  of  the  following  approximations: 

In  ($-1)  =  3-2,  3=2 
or 

In  (3)  =  3-1,  3=1 

For  6  “1.58,  these  are  equally  in  error.  This  maximum  error  is  equal 
to  +.12,  equivalent  to  a  multiplicative  error  in  P  of  1.12. 


We  now  have  two  possible  cases,  specifically. 

In  (g-1)  =:  3-2,  3  > 1.58  ) 

or  > 

In  3  =  3-1,  3  < 1.58  I 

of  course,  equation  (A-9)  implies  a  rough  knowledge  of  the  variable  for 


(A-9) 


which  we  wish  to  solve.  This  is  often  available,  especially  in  a  tabula¬ 
tion  approach,  wherein  P  might  be  held  constant  and  3  calculated  for 
increasing  values  of  N.  Then,  previous  answers  supply  the  necessary 
estimate.  When  no  estimate  is  available,  one  may  choose  at  random  or 
try  both  methods  and  compare.  Depending  upon  the  equation  chosen  for 
(A-9),  one  of  the  following  results. 


Case  1. 


3  >  1.58 


Equation  (A-8)  is  expressed  entirely  in  terms  of  3-  Define  (for  fixed 

P  and  N,  this  is  a  constant) 

ic  N-lnP  +  1/2  InN  -  ln(N-l)  +  (1 -1/2  in  2tt) 

N 


Then,  we  want  to  solve 


Case  2. 


3  -  In  3  =  K 


(A- 10) 


(A-ll) 


3  <  1.58 


Let  3'  =  (3-1).  Define  a  constant 


K'  =  -InP  +  1/2  InN  -  ln(N-l)  -  1/2  In  2ir  =  N(K-1)-1 


(A-12) 


We  then  wish  to  solve 


In  8'  -  8'  -  K' 


In  either  case,  we  have  a  simple  transcendental  equation  which 
does  not  admit  of  direct  solution.  However,  an  iterative  technique  (such 


as  Newton’s)  is  readily  applied. 

Defining  f(S)  -  B-lnS-K, 

f  (3)  “  1-1/8 


(A-14) 


A-4: 


(A-15) 


Newton's  method  for  solving  (A-ll)  is: 

£(3i} 

3i+l  "  3i  "  f'(ei) 

The  initial  estimate  is  obtained  from 
Inf}  =  (8-1)  -  1/2  (g-l)2 

whence 

3o  =  1  +  1/2  (K-l)  (A-16) 

A  similar  method  holds  for  the  solution  of  (A-13) . 

Solution  of  either  (A-ll)  or  (A-13)  yields  an  approximate  inver¬ 
sion  of  the  incomplete  gamma  ratio.  Dependent  upon  the  magnitudes  of 
P  and  N  and  the  desired  accuracy  of  A,  the  result  may  be  sufficient.  If 
the  term  ('vl) ,  neglected  in  the  solution  above,  is  considered  significant, 
the  calculated  value  of  g  may  be  used  as  an  initial  estimate  in  a  more 
exact  iteration  using  equation  (A-8) .  Alternatively,  it  may  be  used  as 
the  initial  estimate  in  the  other  method  we  describe. 

We  present  here  a  few  examples  illustrating  the  accuracy  of 
the  preceeding  method.  The  validity  of  the  results  was  tested  by  comparison 
against  a  table  of  solutions  to  the  inverse  problem  given  in  Ref. 

In  general,  the  larger  N  or  the  smaller  P,  the  more  accurate  the  approxima¬ 
tion. 

1.  from  the  tables,  N  =  100,  A  =  187.248  =>P  =  10-12. 

Using  the  direct  solution,  as  in  equation  (A-7) ,  yields 
P  -  1.005  x  10-12. 

2.  Inverting  the  same  problem,  using  equations  (A-10) ,  (A-ll), 
and  (A-14)  thru  (A-16),  yields  A  ■  187.272.  The  error 
either  here  or  above  is  presumably  insignificant. 
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.  A  similar  inversion,  but  for  smaller  numbers.  From  the 
tables,  N  =  32,  P  =  10“6->A  =  66.3937.  The  calculated 
root  is  A  ■  66.4513. 

4.  Letting  N  get  smaller,  and  P  larger,  we  begin  to  notice 
errors  which  can  be  deemed  significant. 

a)  N  =  10,  P  =  10  ^  —>  A  =  14.206.  The  calculated  value 
is  A  =  14.490. 

b)  Using  N  =  10  and  the  calculated  A,  the  direct  evaluation 
yields  P  =  .83  x  10  \  Using  A  from  the  tables, 

P  -  .92  x  10"1. 

These  examples  point  out  an  interesting  behavior.  Although 
the  inverse  solution  is  derived  from  the  direct  one  by  further  approxima¬ 
tions  and  simplifications,  it  produces  more  accurate  results.  This  is 
simply  a  matter  of  relative  sensitivities.  By  examining  the  tables  in 
Ref.  5.,  a  small  change  in  A  (a  few  percent)  is  seen  to  be  sufficient 
to  change  P  by  an  order  of  magnitude. 

5.  The  last  example  uses  the  smallest  possible  N  and  largest 
P  and  should  represent  worst  performance  for  the  approximation.  Of 
course,  in  this  case,  the  problem  could  be  solved  without  recourse  to 
approximation  as  explained  in  Section  4.2.2.  Nonetheless,  let  N  *  2  and 
P  ■  10  \  The  true  solution  is  A  =  3.8897;  the  calculated  inverse  is 

A  *  3.6636.  Evaluating  directly  with  the  true  root  yields  P  =  1.58  x  10 

We  discuss  now  an  alternative  computational  method  for  evaluating 
and/or  inverting  the  incomplete  gamma  ratio.  This  is  the  method  described 
by  Pachares  in  Ref.  5.  We  use  it  here  primarily  to  calculate  detection 
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probabilities,  where  A  <  (N-l) ,  and  the  previous  approximation  breaks 
down.  Pachares  used  it  to  Invert  the  false-alarm  equation  to  calculate 
the  values  of  A  in  his  tables.  To  ensure  accuracy  in  determination  of 
the  performance  results  given  in  Section  4.2.3,  we  also  do  this.  How¬ 
ever,  the  value  of  A  given  by  the  preceeding  method  is  used  as  an  initial 
estimate  in  Pachares'  iterative  approach.  This  is  more  convenient  to  use 
than  the  initial  estimate  he  considers. 

As  discussed  in  Section  4.2.2,  P(N,A)  may  be  evaluated  exactly 
as  a  sum  of  N  terms. 


P(N,A) 


k=0 


This  may  be  re-structured  to  simplify  the  computations. 


Define 


*<"■«  ■  M5-  2  fr 


k=0 


(A-17) 


Then, 


-A  .N-l 

P(N,A)  -  e/M_T—  S(N,A) 


(N-l) I 

By  suitable  manipulation,  S(N,A)  can  be  expressed  in  a 
recurrence  relationship , 

S(N,A)  =  1  +  S(N-1,A) 

S(1,A)  =  1 

Although  tedious,  evaluation  of  this  involves  terms  whose  magnitudes  are 

easy  to  handle.  S  does  not  get  to  be  very  large - as  an  example,  N  *  100 

and  A  -  187.248  (==>P  =  10-^)  yield  S  =  2.096.  Then,  P  can  be  rather 
readily  evaluated  from  (A- 18) .  The  only  approximation  that  need  be  made 
is  Stirling's  for  the  factorial.  Using  the  form 


(A- 18) 


(A-19) 


In  <N.')  =s  (N+l/2)  InN  -  N  +  1/2  In  (2n)  + 


12N 


results  in  an  accuracy  of  at  least  9  places  for  N  >7 

A -7 


1 

360N3 


To  invert,  we  wish  to  calculate  the  root  of 
f (A)  A  p(N,A)  -  P  -  0 

Newton’s  method  for  the  solution  of  this  involves  iterations  of  the  form 


i+1 

Since  (by  original  definition) 

P(N 


A  f(Ai> 

i  "  FTO 


r  e“t  t1*-1 
♦A>  =  I  fefr-  dt 


we  have 


d  e-A  aN-1 

3a  -  - 


(N-l) ! 


The  iteration  then  reduces  to 


Ai+1  “  Ai  +  S<N’V  “  A(N»Ai> 


where 


A(N,  A)  = 


(N-l) ! 


e“A  A11'1 


(A-20) 

(A-21) 


A  (N,A)  is  readily  calculated  and  is  of  the  same  order  of  magnitude  as 
S.  Thus,  evaluation  of  equation  (A-20)  is  a  practical  technique  for 
inversion.  Using  a  good  starting  value,  such  as  given  by  the  solution 
of  (A-ll)  or  (A-13) ,  only  a  few  iterations  are  required  to  reach  a 


solution. 
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ABSTRACT 

This  report  is  a  logical  companion  to  Technical  Note  76-01, 
"Detection  of  a  Message  Redundantly  Transmitted  Over  a  Fading  Channel". 
In  this  report,  we  analyze  the  performance  of  a  would-be  interceptor, 
comparing  his  ability  to  detect  the  message  versus  the  capability  of 
the  intended  receiver,  using  a  measure  termed  "recognition  differential" 
Performance  curves  for  two  interceptors,  which  form  upper  and  lower 
bounds  to  the  performance  of  any  interceptor,  are  presented. 
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SECTION  1 

INTRODUCTION  AND  SUMMARY 

In  a  companion  report  (Ref.  1)  we  presented  a  communications 
problem  and  an  approach  to  its  solution.  Briefly,  the  requirement  was  to 
transmit  with  high  reliability  the  logical  equivalent  of  one  bit  of 
information.  The  constraints  included  a  lack  of  synchronization  between 
transmitter  and  receiver,  the  use  of  a  fading  dispersive  channel  and  the 
need  to  be  covert.  A  message  structure  and  optimum  receiver  processing 
were  proposed.  Idealized  performance  results  (expressed  in  terms  of 
detection  probability  versus  signal-to-noise  ratio)  were  derived. 
Modifications  to  these  results  to  account  for  the  effects  of  the  dispersive 
nature  of  the  channel  were  also  discussed. 

In  this  paper,  we  analyze  the  rest  of  the  problem---i.e. ,  the 
covertness  of  the  proposed  signalling  scheme.  Two  potential  interceptors 
are  described  and  their  performances  are  derived  in  a  manner  analogous  to 
that  for  the  intended  receiver.  These  levels  of  performance  are  then 
contrasted  against  that  of  the  intended  receiver.  The  difference  is 
termed  recognition  differential  and  is  a  measure  of  the  covertness  of 
the  chosen  message  structure. 

As  explained  in  our  earlier  report,  the  channel  imposed 
distortions  are  such  that  incoherent  processing  is  all  that  is  available 
to  the  intended  receiver.  This  forces  the  use  of  techniques  similar  to 
those  which  would  be  employed  by  an  unauthorized  receiver,  and  much  of 
the  processing  gain  available  in  more  conventional  LPI  systems  is  lost. 
Here,  we  show  that  relatively  long  (e.g.,  100  tones)  messages  are  needed 
to  achieve  a  reasonable  (6  to  10  dB)  value  of  recognition  differential. 
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SECTION  2 


STATEMENT  OF  THE  PROBLEM 


2.1  THE  MESSAGE 

In  this  section  we  summarize  the  format  of  the  transmitted 
message  as  detailed  in  Ref.  1.  The  message  is  assumed  to  consist  of 
a  total  of  N  tones,  where  N  is  a  design  parameter.  Each  tone  has  the 
same  basic  pulse  shape  and  is  centered  in  a  unit  area  of  duration  T  and 
width  of  Af.  Time  is  correspondingly  quantized  into  intervals  of  length 
T,  frequency  into  intervals  of  width  Af ,  forming  a  time -frequency  matrix. 
The  message  tones  are  located  in  the  matrix  according  to  a  frequency¬ 
hopping  pattern. 

Specifically,  we  consider  a  message  of  length  N  located  within 
an  NxN  matrix  such  that  at  any  time  only  one  frequency  is  used,  and 
each  frequency  is  used  at  only  one  time.  The  particular  choice  of  tone 

i 

locations  used  in  the  matrix  is  controlled  by  the  message  pattern  or 
key.  This  is  assumed  known  to  the  intended  receiver  but  not  to  an 
unauthorized  or  intercept  receiver.  It  is  this  knowledge  which  enables 
the  intended  receiver  to  detect  a  message  at  a  lower  received  signal -to- 
noise  level  than  the  interceptor  can,  and  it  is  this  "processing  gain" 
which  (as  a  function  of  the  message  length  N)  we  wish  to  examine. 
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2.2  ASSUMPTIONS 


Several  assumptions  regarding  the  processing  done  by  a  would-be 
interceptor  are  briefly  stated  in  this  section.  In  general,  these 
parallel  those  made  concerning  the  authorized  receiver  in  our  previous 
analysis.  Following  this,  we  describe  the  general  form  that  our  analyses 
will  take. 

Due  to  the  Rayleigh  fading  channel,  the  received  signal  tones 
are  samples  of  a  complex  Gaussian  process.  It  is  assumed  that  the  time 
and/or  frequency  separation  between  active  tones  in  the  message  pattern 
is  sufficient  to  ensure  that  they  are  all  mutually  independent.  It  is 
further  assumed  that  they  are  all  received  at  the  same  energy  level. 

The  background  noise  is  taken  to  be  additive  Gaussian,  independent  and 

identically  distributed  across  the  time -frequency  matrix. 

2  2 

The  noise  variance  is  <r  ,  that  of  the  signal  is  o-  ;  their 

s 

ratio  is  the  signal -to -noise  ratio  denoted  7,  which  is  assumed  constant 
for  all  of  the  received  message  tones.  Letting  u  be  the  response  out 
of  an  arbitrary  time-frequency  resolution  cell,  we  then  have  two  basic 
sub -hypotheses : 

hQ  :  u  e  complex  N(0,cr^)  ,  noise  only  J 

h^  :  u  6  complex  N(0,o^/a)  ,  signal  +  noise  \  (2-1) 

with  a  ^  y+y  \ 

Although  the  above  may  seem  an  obvious  starting  point,  it 
tacitly  makes  several  "worst-case"  assumptions  regards  the  interceptor. 
First,  it  assumes  that  he  is  forming  the  proper  time-frequency  matrix 
at  his  receiver.  That  is  he  knows  the  pulse  length  and  width  and 
processes  his  input  signal  accordingly.  Further,  he  is  synchronized 
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in  both  time  and  frequency  with  the  message  and  there  is  no  run-out  due 
to  Doppler  shift.  We  will  also  assume  that  the  interceptor  knows  the 
total  number  of  active  tones  (N)  and  the  dimensions  of  the  matrix 
containing  them. 

The  assumption  of  synchronization  is  one  that  was  also  made 
in  our  earlier  analysis.  There,  we  described  how  to  model  the  more 
realistic  unsynchronized  situation  by  decreasing  the  effective  SNR  per 
tone.  That  discussion  applies  equally  well  to  the  intercept  receiver, 
in  fact,  when  comparing  SNRs  for  the  two  receivers,  this  extra  factor 
washes  out. 

Whether  or  not  an  interceptor  would  in  reality  know  the  form 
of  the  message  (T;  Af;  N)  is  not  within  the  scope  of  this  study.  We 
choose  these  assumptions  using  the  standard  rationalization  behind  any 
worst-case  analysis.  However,  in  interpreting  the  results,  their  basic 
pessimistic  nature  should  not  be  forgotten. 

For  simplicity  of  description,  we  will  picture  the  interceptor 
as  forming  and  storing*  an  NxN  array  of  responses  exactly  as  the  authorized 
receiver  would.  Incoherent  (square-law)  processing  will  always  be  used, 
so  only  the  magnitude-squared  of  each  resolution  cell  need  be  saved. 

Only  two  cases  will  be  of  interest.  The  noise-only  hypothesis  (Hq) 
assumes  that  none  of  the  N  tones  in  the  message  are  within  the  matrix 
currently  being  examined;  the  signal  hypothesis  (H^)  assumes  that  they 
all  are. 


As  explained  in  Section  2.4,  this  is  not  really  necessary. 
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According  to  some  algorithm,  the  interceptor  will  process 
the  responses  in  his  received  time-frequency  matrix  and  form  therefrom 
a  decision  variable.  Just  as  for  the  authorized  receiver,  this  will  be 
compared  against  a  threshold  to  determine  whether  a  message  is  present 
or  not.  The  setting  of  this  threshold  will  be  chosen  to  yield  the 
desired  per-decision  false-alarm  probability  (P  .  ).  The  signal -to -noise 
ratio  required  by  the  interceptor  to  achieve  probability  of  detection  =  0.5 
is  computed.  This  value  is  then  compared  to  that  required  by  the  authorized 
receiver  to  similarly  have  a  fifty  percent  detection  probability.  The 
difference  (in  dB)  between  these  signal -to-noise  ratios  is  termed 
"recognition  differential"  and  is  a  direct  measure  of  the  processing 
gain  that  the  intended  receiver  has  over  an  interceptor. 

In  general,  an  interceptor's  received  SNR  will  differ  from 
that  of  the  authorized  receiver.  Three  possible  causes  are 

1)  Different  range  from  the  transmitter. 

2)  Different  bearing  relative  to  the  transmitted  beam 
center  (MRA-Main  Response  Axis) . 

3)  Different  directionality  index  of  the  receiving  array. 

Any  of  these  can  result  in  a  received  SNR  greater  or  less  than  that  of 
the  intended  receiver.  Here,  we  do  not  concern  ourselves  with  particular 

combinations  of  these  three  possibilities - we  merely  derive  the  required 

net  effect.  Given  a  particular  operational  scenario,  recognition  differ¬ 
ential  may  then  be  expressed  in  terms  of  one  of  these  parameters.  As 

an  example,  if  an  interceptor  were  presumed  to  use  a  receiving  array  of 
the  same  directivity  and  was  at  the  same  angle  relative  to  MRA  as  the 
intended  receiver,  then  the  recognition  differential  could  be  equated 
to  a  "range  advantage"  for  the  intended  receiver. 
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2.3  POSSIBLE  INTERCEPTORS 

For  the  intended  receiver,  it  was  obvious  that  only  those  cells 
in  the  time -frequency  matrix  (as  determined  by  the  key)  corresponding 
to  signals  should  be  considered  in  the  formation  of  a  decision  as  to  the 
presence  of  a  message.  The  interceptor  does  not  have  such  an  obvious 
strategy,  in  fact,  several  forms  of  processing  could  equally  well  be 
postulated,  each  representing  differing  degrees  of  sophistication  and/or 
knowledge  of  the  message  structure.  We  briefly  discuss  some  of  these 
and  then  select  two  for  analysis. 

The  most  simplistic  form  of  intercept  analysis  examines  only 
one  time  slice  or  frequency  cell.  As  discussed  in  Section  2  of  our 
earlier  report,  this  leads  to  a  message  structure  of  at  most  one  tone 
per  time  or  frequency  in  order  to  minimize  the  performance  for  this 
interceptor.  Although  an  interceptor  may  also  use  more  complex  strategies, 
we  take  this  as  a  baseline  system  against  which  covertness  must  be 
considered.  Hence,  we  choose  our  basic  message  as  indicated  to  optimize 
performance  against  this  basic  interceptor.  In  order  to  evaluate  more 
significant  threats,  we  now  consider  ways  in  which  an  interceptor  can 
process  over  the  entire  NxN  matrix  in  an  attempt  to  detect  the  presence 
of  a  message. 

We  will  choose  two  interceptors  which  may  be  construed  as 

bounding  the  performance  of  any  sophisticated  yet  realistic  threat.  The 

first  is  a  traditional  energy  detector  (denoted  ED) .  This  interceptor 

* 
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8 imply  sums  all  of  the  N  responses  in  the  matrix.  This  decision  variable 
is  then  compared  to  a  threshold.  Although  the  ED  is  a  sophisticated  processor 
in  that  it  makes  a  decision  after  examining  the  entire  received  NxN 
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matrix,  it  represents  the  most  elementary  form  of  processing  such  a  matrix. 
Thus,  we  view  its  performance  as  being  a  lower  bound  to  the  performance 
of  an  arbitrary  interceptor. 

We  now  additionally  assume  that  the  interceptor  knows  that  the 
message  is  structured  so  that  only  one  tone  is  transmitted  at  any  time. 

With  this  knowledge,  a  better  form  of  processing  is  possible.  Specifically, 
for  each  time  slice,  the  greatest  received  power  level  in  the  N  frequency 
cells  is  selected  and  saved.  The  sum  of  these,  accumulated  over  N  time 
slices,  forms  the  decision  variable.  We  refer  to  this  as  the  "greatest 
of"  receiver,  denoted  GOF. 

Knowing  that  only  one  tone  is  transmitted  at  a  time  defines 
possible  patterns.  Further  specifying  (which  we  do  not  consider  here) 
that  no  frequencies  are  repeated  lowers  this  number  to  Nl.  An  optimum 
interceptor  would  form  a  decision  variable  corresponding  to  each  of 
these  possible  message  patterns.  Without  possession  of  the  key,  the 
interceptor  must  consider  each  of  these  as  representing  a  possible 
message.  The  most  likely  message  is  simply  the  one  with  the  greatest 
associated  decision  variable.  An  optimum  interceptor  would  test  this 
most  likely  message  to  decide  if  any  message  were  present.  But,  by 
maximizing  the  likelihood  ratio  in  each  time  slice,  the  GOF  directly 
extracts  this  maximum  likelihood  term.  Thus,  the  GOF  is  equivalent 
to  the  optimum  form  of  processing,  and  is  clearly  within  the  con¬ 
straints  of  practicality.  We  consider  the  GOF  as  representing  the 
optimum  interceptor,  and  we  interpret  its  performance  as  being  an 
upper  bound  to  that  of  any.  interceptor. 
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Having  concocted  an  interceptor  optimized  against  our  message 
structure,  it  is  tempting  to  consider  changing  the  message  structure. 

For  instance,  if  in  one  time  slice  two  tones  were  transmitted  and  then 
none  in  the  next,  the  GOF  would  clearly  suffer  a  decrease  in  performance. 
However,  such  reasoning  is  circular  as,  having  postulated  this  new 
message,  one  could  in  turn  define  a  new  intercept  strategy,  optimized 
against  it.  Thus,  we  fix  our  message  structure  as  described  and  merely 
analyze  how  well  these  two  interceptors  perform  against  it. 

We  might  further  consider  other  intercept  strategies.  One 
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such  might  be  termed  "threshold  counting."  Here,  each  of  the  N  responses 
is  compared  against  a  fixed  threshold  set  by  the  noise  variance  and 
the  number  of  exceedances  is  counted.  This  total  is  then  compared  against 
another  threshold  and  a  decision  as  to  the  presence  of  a  message  is 
made.  This  technique  is  related  to  the  cognitive  processes  of  a  trained 
operator  who  is  viewing  the  responses  on  a  screen  and  looks  for  an 
unusually  large  number  of  "blips."  It  is  probably  less  useful  as  a 
model  of  a  technique  which  would  be  implemented  purely  electronically. 
However,  our  subsequent  analysis  shows  that  our  upper  and  lower  bounds 
to  the  performance  of  an  arbitrary  interceptor  are  relatively  close. 

Thus,  we  do  not  analyze  in  detail  the  threshold  counter  or  any  other 
possible  intercept  scheme. 
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PROCESSING  CONSIDERATIONS 


In  this  section  we  discuss  briefly  some  significant  differences 
between  an  interceptor  and  an  authorized  receiver.  The  authorized 
receiver  has  an  advantage  in  that  he  alone  possesses  the  message  key. 
However,  this  also  penalizes  him.  As  explained  in  Ref.  1,  he  must 
examine  the  received  time -frequency  matrix  every  time  slice  in  order  to 
detect  the  message.  Further,  multiple  Doppler  hypotheses  are  required 
to  properly  align  the  received  tones  with  the  pattern. 

Since  the  interceptor  does  not  have  the  pattern,  he  does  not 
care  about  the  fine  alignment  of  his  time-frequency  matrix  with  the 
message.  Not  only  need  he  not  form  multiple  Doppler  hypotheses,  he 
need  not  examine  every  matrix.  For  instance,  examination  of  only  every 
tenth  matrix  results  in  missing  at  most  5  tones.  For  a  message  consisting 
of  128  tones,  this  represents  a  loss  in  signal  energy  of  only  .2  dB, 
which  can  be  considered  as  being  negligible.  Thus,  the  interceptor  need 
make  decisions  much  less  frequently  than  the  intended  receiver.  We  also 
note  that,  for  both  the  ED  and  GOF,  successive  decision  variables  will 
be  highly  correlated. 

The  significance  of  all  this  is  that,  in  order  to  have  a  false- 
alarm  interval  the  same  as  that  for  the  intended  receiver  (which  seems 
the  most  reasonable  manner  in  which  to  compare  their  performance),  the 
interceptor  can  operate  at  a  higher  false-alarm  probability.  This  is 
equivalent  to  a  lowering  of  his  threshold,  thereby  increasing  his  detection 
probability.  The  difference  between  operating  P^s  can  be  one  or  two 

orders  of  magnitude.  In  our  comparisons,  we  will  assume  that  the  intended 
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receiver  is  operating  at  ■  10  and  the  interceptor  at  10  . 

\ 
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SECTION  3 


OUTLINE  OF  ANALYSES 

Our  analysis  of  an  intercept  receiver  is  somewhat  different  in 
approach  from  our  earlier  analysis  of  the  performance  of  the  intended 
receiver.  It  involves  initially  making  two  assumptions  (valid  for  large 
N)  and  then  determining  correction  terms  for  smaller  N.  As  this  approach 
is  common  to  our  discussion  of  both  the  Energy  Detector  and  the  Greatest 
of  interceptors,  we  outline  it  here. 

We  consider  two  idealized  hypotheses,  and  H^,  corresponding 

respectively  to  receiving  noise  only  and  noise  plus  all  N  message  tones. 

The  calculated  decision  variable  is,  in  general,  denoted  as  v.  To 

emphasize  the  dependence  of  the  statistics  of  v  upon  the  particular 

hypothesis,  we  will  write  v^  with  j  equal  to  either. 0  or  1.  Similarly, 

where  it  is  necessary  to  explicitly  state  the  number  of  terms  used  in 

the  combination  forming  v  (i.e.,  the  order  of  the  random  variable)  this 
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shall  be  represented  as  an  argument,  e.g.,  v^  (N  ) . 

We  note  that  for  the  ED  and  GOF,  v^  is  formed  as  the  sum  of 
many  independent  and  identically  distributed  (lid)  random  variables.  We 
thus  approximate  v^  by  a  normal  random  variable  with  appropriate  mean 
and  variance.  This  allows  us  to  calculate  a  decision  threshold  which 
closely  approximates  the  ideal  threshold  corresponding  to  a  specified 
as  follows: 

2 

First,  the  mean  (u)  and  variance  (o^  )  of  Vq  are  calculated. 
From  a  tabulation  of  the  standard  (i.e.,  zero  mean  and  unit  variance) 
normal  probability  function,  such  as  in  Table  3.1,  the  value  \ 
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corresponding  to  is  found.  Then,  the  approximate  threshold  for  the 
interceptor  is  obtained  as 

A  =  P  +  *0q  (3-1) 

With  the  threshold  established  as  just  described  we  next 

consider  the  signal  present  hypothesis  .  We  determine  the  required 

SNR  to  achieve  an  intercept  probability  of  0.5.  This  point  corresponds 

to  the  SNR  for  which  the  median  of  v^  equals  the  threshold.  We  make 

our  second  approximation,  which  is  that  the  mean  of  and  the  median 

of  Vj  coincide.  This  relationship  is  strictly  valid  for  symmetrical 

distributions  and  is  closely  approximated  for  the  high -ordered  distributions 

of  interest  here.  To  have  a  fifty  percent  detection  probability  is 

thus  equivalent  to  specifying  that 

E  [vj]  *  A  (3-2) 

The  expected  value  of  v^  may  be  expressed  as  a  function  of 

* 

the  signal-to-noise  ratio.  This  allows  us  to  invert  Equation  3-2  to 
solve  for  SNR  as  a  function  of  A.  Using  Equation  3-1,  we  then  eliminate 
A,  resulting  in 

E(vl3  =  EIV0]  +  A\/varfvo]  (3"3) 

This  is  solved  for  SNR  as  a  function  of  N,  completing  our  initial 

i 

analysis . 

We  wish  to  produce,  as  our  measure  of  performance,  a  curve 
showing  required  SNR  vs  message  length  N.  The  results  determined  from 
Equation  3-3  are  onlv  asymptotically  correct  for  large  N.  By  more 
precisely  calculating  points  for  smaller  N,  we  can  modify  this  curve  to 
eliminate  the  errors  introduced  by  our  assumptions.  By  their  definition, 
both  the  ED  and  the  GOF  are  equivalent  to  the  intended  receiver  for  the 

I 


I 
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degenerate  case  of  N  = 1.  Thus,  the  exact  performance  is  already  known 
for  N  =  l.  More  accurate  calculations  for  a  few  other  values  of  N  then 
allows  a  smooth  curve  to  be  fitted  to  the  asymptotic  result. 

Of  our  two  assumptions,  the  first  represents  the  greatest  source 
of  error.  Here,  the  central  limit  theorem  is  used  to  approximate  the 
distribution  of  v^  in  its  tails,  and  it  is  here  that  the  limiting  form  is 
a  poor  approximation  for  small  N.  More  detailed  discussions  of  accuracy 
and  corrections  are  presented  in  the  particular  analysis  following. 

The  second  assumption  is  not  regarded  as  introducing  as 
significant  an  error.  For  one  reason,  it  is  a  weaker  assumption  than 
the  first.  The  condition  that  the  mean  and  median  coincide  is  satisfied 
by  any  random  variable  with  density  symmetric  about  its  mean,  whether 
it  be  normal  or  not.  Further,  the  exact  shape  of  the  tails  of  the 
density  of  v^  have  little  effect  on  the  desired  condition.  As  an  example, 
for  a  chi-squared  density  of  order  greater  than  21,  the  difference 
between  the  mean  and  the  median  is  less  than  .1  standard  deviation. 


Also,  some  non-symmetric  examples  can  be  concocted. 
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SECTION  4 


THE  ENERGY  DETECTOR 


The  energy  detector  forms  a  decision  variable  by  adding  the 

magnitude-squared  of  all  the  responses  in  the  matrix.  This  is  normalized 
2  2 

by  the  factor  2o  ,  where  a  is  the  (presumed  known)  variance  of  the  noise 
process  and  then  compared  against  a  fixed  threshold  A.  The  decision 
variable  is  in  general  denoted  by  v,  more  specifically,  by 


vj  «>  -  rj  S  h 


(4-1) 


where  the  subscript  j  specifies  the  hypothesized  distribution  of  the 
responses  ju^J.  and  the  argument  R  is  the  number  of  terms  used  in  forming 


For  an  NxN  matrix,  the  two  hypotheses  are  given  by 
Hq!  all  N2  terms  are  purely  noise 

(N2-N)  terms  are  noise;  N  are  signal  plus  noise 


(4-2) 


2 

Under  H^,  the  ju^J  are  i.i.d.  and  v^(N  )  has  a  standard  Gamma  distribution 


of  order  N  .  We  thus  have 


E  O0(N2)]  -  N 
Var  [vQ(N2)]  =  N2 


(4-3) 


Under  H^,  the  summation  indicated  in  Equation  (4-1)  may  be 


re-written  as 


ri<1,2)  =  r?  E  |ui|2  +  71  E  |ui| 

2  0  *  1  2a  i=l  *  1 


where  the  sub-hypotheses  h^  and  h^  determining  the  distribution  of  the 
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individual  u^  are  as  defined  in  Equation  (2-1) .  We  then  have  two 

independent  summations.  The  first  is  simply  a  standard  Gamma  variate 
2 

of  order  (N  -N) .  The  second  is  also  a  Gamma  variate,  or  order  N,  but 

is  non-standard— -i.e. ,  it  has  a  shape  parameter  not  equal  to  unity. 

2 

Although  v^(N  )  is  the  sum  of  two  Gamma  variates,  since  they  have 

differing  shape  parameters,  the  reproductive  property  of  the  Gamma 

distribution  does  not  hold.  Thus,  their  sum  is  not  a  Gamma  variate  and 

its  density  must  be  determined  by  convolution. 

By  normalizing  the  second  sum  by  the  factor  1/a ,  it  can  be 

2 

written  in  terms  of  a  standard  Gamma  variate.  Then,  v^(N  )  may  be 
expressed  as  a  function  of  two  variables  formed  under  Hq, 

Vl(N2)  =  v0(N2-N)  +  i  vq(N)  (4-4) 

From  this,  we  directly  have  the  moments  of  v1 : 

E  [vx(N2)] 

and , 

var  \y1  (N2)] 

or, 

E  [y1  (N2)] 
var  [vx  (N2)] 

Inserting  (4-7)  and  (4-3)  into  Equation  (3-3),  we  find  the 
SNR  required  for  the  ED  to  achieve  P^= . 5  to  simply  be  given  by 

Y  ■=  X  (4-9) 


=  E  [v0(N2-N)]  +  i  E  [v0(N)] 


=  (N  -N)  +  (1-rY)  (N) 

(4-5) 

=  N2  +  YN 

-  (N2-N)  +  (1+Y) 2  (N) 

(4-6) 

=  N2  +  2YN  +  Y2N 

=  E  [v0(N2)]  +  YN 

(4-7) 

A  2 

=  var  [vq(N  )]  +  2YN  +  Y  N 

(4-8) 
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This  rather  surprising  result  states  that,  as  long  as  N  is  large  enough 
for  our  approximations  to  be  valid,  Y  is  independent  of  N. 

As  explained  in  Section  2.4,  we  shall  specifically  consider 
an  interceptor  operating  at  a  false-alarm  level  of  10  We  then  have, 
from  Table  3-1, 

Y  =  4.265  or  Y=  6.30  dB  (4-10) 

This  is  shown  as  the  solid  line  in  Figure  4-1. 

This  result  is  asymptotically  correct  for  large  N.  We  now 

determine  more  exact  results  for  small  N.  First,  we  note  that  the 

primary  source  of  error  is  in  the  setting  of  the  threshold  using  the 

2 

values  in  Equation  (4-3).  For  the  ED,  the  exact  threshold  A(N  )  is 

easily  enough  determined.  Using  the  tables  given  in  Reference  2,  values 

-5 

of  the  threshold  A  yielding  a  of  10  are  found.  The  requisite  SNR 

is  then  calculated  using  Equations  (3-2)  and  (4-5),  viz., 

N2+YN  =  A(N2)  (4-11) 

These  results  are  shown  as  the  dashed  line  in  Figure  4-1. 

Finally,  we  examine  the  error  introduced  by  equating  the 
2 

median  of  v^(N  )  to  its  mean.  This  necessitates  determining  the  density 
of  Vj  by  convolving  the  densities  of  the  random  variables  indicated  in 
Equation  (4-4).  This  is  then  integrated  to  determine  the  fiftieth 
percentile  point.  Equating  this  to  the  exact  threshold  will  eliminate 
all  approximations  and  the  determined  answer  will  be  an  exact  level  of 
performance  for  the  ED. 

For  N=l,  the  exact  SNR  for  PD=.5  is  the  same  as  for  the  intended 
receiver  and  is  readily  obtained  from  our  earlier  analysis. 


ia 


FIGURE  4-1  PERFORMANCE  OF  THE  ED:  SNR  FOR  PD  -  .5  @  PFA 


From  Figure  4-3  of  Reference  1,  we  have  Y»11.9  dB.  For  N=2,  the  convolution 
and  integration  are  feasible,  although  already  awkward.  The  required  SNR 
is  calculated  to  be  9.50  dB.  Exact  values  could  be  found  for  larger  N, 
but  the  effort  does  not  appear  justified,  as  the  corrections  are  small. 

A  sufficiently  smooth  curve  can  be  fitted  through  these  two  exact  results 
to  the  asymptotic  result.  This  is  the  dashed  and  dotted  curve  in 


Figure  4-1 
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SECTION  5 

THE  GREATEST  OF  RECEIVER 


Out  of  each  time  slice,  the  "greatest-of"  receiver  extracts 

the  greatest  (magnitude -squared)  response.  This  is  done  for  the 

appropriate  number  of  time  slices.  The  selected  terms  are  summed  and 

2 

normalized  by  the  factor  2 o  to  form  the  decision  statistic.  This  is 
then  compared  against  a  fixed  threshold  A. 


As  we  did  before,  we  will  denote  the  decision  statistic  as  v, 

or,  more  specifically  as  Vj (K) ,  where  j  specifies  the  hypothesis  and  K, 

the  number  of  terms  summed  to  form  v.  Whereas  for  the  ED  v  was  formed 

2 

by  directly  adding  the  j |u^|  } >  we  now  £ orm  v  as  the  sum  of  random 
variables  which  are  functions  of  the  {  These  functions  represent 

the  "greatest  of"  operation.  We  denote  the  individual  terms  by  w,  or 
more  specifically,  by 


w.(K')  &  max 

3  2o* 

where  the  subscript  j  specifies  the  hypothesis  controlling  the  distribution 
of  the  {uj  and  the  argument  K'  is  the  number  of  responses  from  which 
the  greatest  is  extracted.  Note  that  for  convenience  in  the  analysis 
the  normalizing  factor  2o  (where  q  is  the  known  noise  variance)  is 
applied  directly  to  the  {  w  }.  The  decision  statistic  is  then  defined  as 

K  . 

Vj  (K)  £  (wj(K,>)i  <5“2> 

The  arguments  K  and  K'  correspond  respectively  to  the  number 
of  time  slices  and  the  number  of  frequency  cells  in  the  time/frequency 
matrix.  .  s  described  in  Section  2.2,  we  are  considering  the  special 
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case  of  an  NxN  matrix,  hence,  K  and  K1  are  equal.  In  a  more  general 
formulation,  however,  one  would  need  to  properly  distinguish  them.  Here 
we  set  them  both  to  N. 

As  we  are  assuming  the  signal  and  noise  terms  throughout  the 
matrix  to  all  be  mutually  independent,  we  clearly  have  the  {u^}  defining 
w  in  Equation  (5-1)  all  mutually  independent;  similarly  so  are  the  |w} 
defining  v  in  Equation  (5-2).  The  distribution  of  the  individual  |u^| 
in  (5-1)  is  determined  by  the  choice  of  hypothesis: 

Hq!  all  N  of  the  |u^}  are  noise  (hg) 

H^:  (N-l)  of  the  ju^  are  noise  (h^)  ;  one  is  signal  +  noise 

where  the  sub -hypotheses  h^  and  h^  are  defined  in  Equation  (2-1).  Thus, 
within  one  time  slice,  the  ju^,  although  independent,  are  not  necessarily 
identically  distributed.  However,  from  one  time  slice  to  the  next,  the 
individual  greatest  of  terms  appearing  in  the  summation  of  (5-2)  are 
iid. 

We  first  determine  the  probability  density  function  of  w  under 
Hq  and  Hj.  This  is  done  by  finding  the  distribution  function  and 
differentiating. 

Examining  an  individual  response  under  the  two  sub -hypo theses 
of  Equation  (2-1) ,  we  calculate  the  probability  that  this  response  is 
less  than  some  value  w.  The  normalized  magnitude  squared  is  simply  an 
exponentially  distributed  random  variable,  hence 


By  the  assumed  independence  of  the  {uj  ,  we  have 


max 


n2H  ■  fr?^i“i|2  <v* 

H, 


(5-5) 


H, 


This  specifies  the  distribution  function  of  w^(N),  which  we  denote  as 
Fj^(w).  Substitution  into  Equation  (5-5)  from  Equations  (5-4)  and  (5-3) 
then  yields: 


W")  -  (1-e">N 

„  ,  .  ,,  -w.N-1  ..  -aw. 

f1>n(w)  •=  (1  -  e  )  (1-e  ) 


(5-6) 

(5-7) 


The  density  functions,  denoted  by  f  rather  than  F  are  readily  found  to 
be: 

_ _ -  XT  1  _ _ 

(5-8) 


*  ✓  \  -w.N-1 .  -w. 

f0jN(w)  -  N(1  -e  )  (e  ) 


.  ,  .  ....  -w.N-2  -w.,  -aw.  ,  -aw,,  -w.N-1  . 

fl>N(w)  =  (N  -  1)  (1  -  e  )  e  (1-e  )  +  ae  (1-e  )  (5-9) 

The  density  of  w  under  may  alternatively  be  expressed  in  terms  of 
the  density  under  HQ: 


fl,N(w)  =  f0,N-l(w)  "  (N-l)  (1  -  e-w)N"2e“(We'w 
+  a(l— e‘w)N‘1e-°w 


(5-10) 


Next,  we  evaluate  the  moments  of  w.  Details  are  presented  in 
Appendix  A.  We  first  show  that 

-  2E  Q  <-»k'1  (c)2 

E[»l(K)]  -  E[.0(H-»]  +  £.  Q  (-«k  (^) 

e[“i2<n>]  ■  e[»0Vd]  +2^  (N‘1)  <-«k  (k+s)2 
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Simplification  of  these  leads  to  the  following  results.  The  mean  and 


variance  of  w^  are  unexpectedly  simple  and  elegant  ih  form. 


E[w0oo]  -  E1  (i) 


(5-11) 


[»0W)]  -  £  (i) 


E[Wl(N)-]  -  E  [w0(N-D]  +  B  (at,N) 


(5-12) 


(5-13) 


It  does  not  appear  possible  to  simplify  the  expression  for  the  variance 
of  w^,  however,  we  do  not  really  need  to  know  this. 

The  expression  for  the  mean  of  Wq  is  the  well  known  harmonic 
series,  which  is  easily  shown  to  be  divergent.  Thus, 

E[»0<S)]  -* 

The  variance  is  a  lesser  known  series;  it  converges  and  can  be  succinctly 
expressed  in  the  limit,  as 

[w„w]  -  £ 

The  function  B(a,N)  appearing  in  the  mean  of  w.  is  the  Beta  function. 


defined  by. 


l(ct,N)  =  B(N,a)  =  J(1  -  t)N_1  t*"1  dt 
0 


This  may  be  expressed  in  terms  of  Gamma  functions, 


B(a  jj)  . 

V  '  r(a  +iN) 

which,  for  computational  purposes,  is  written  as 


(5-14) 


m  (N  -  1)  (N  -2) 

B(o,N)  -  (N.1+a)  x  (N.2+a)X 


(1  +  a) 


a  <5-15) 


Values  of  B(a,N),  presented  against  Yin  dB  (from  whence  a  is  determined), 
are  shown  in  Table  5-1. 


5-4 


jam 


ro 

TT 

r- 

c 

-4 

o 

x 

* 

•ft 

P- 

•tr 

CN 

o 

-f  a- 

r- 

o 

''O 

•H 

in 

*n 

CN 

JO 

n 

X 

X 

JO 

■o 

JO 

X* 

*2 

r* 

•4“ 

CN 

rH 

'O 

N 

Tv 

JO 

rH 

H 

CN 

*N 

CO 

4* 

m 

r- 

X 

JO 

CO 

JO 

CO 

r- 

CN 

Tv 

CO 

CN 

o 

J0 

•■H 

X 

o 

TV 

CN 

X 

X 

'N 

JO 

CN 

Tv 

n 

CN 

X 

CN 

o 

O 

o 

c 

O 

o 

x 

o 

o 

■h 

•H 

CN 

ro 

JO 

r- 

rH 

r- 

rr 

JO 

Tv 

r* 

H 

CN 

CN 

CN 

X 

r* 

Tv 

CT 

r- 

in 

— 

«H 

CC 

CO 

r- 

JV 

JO 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

H 

CN 

3 

4 

X 

T> 

CN 

X 

rH 

*n 

JO 

X 

JO 

JO 

rH 

4* 

X 

X 

rH 

CO 

X 

r 

Tv 

JO 

rH 

rH 

CN 

CN 

ro 

4 

in 

O' 

CN 

JO 

Tv 

*4 

*H 

Tv 

cr 

CN 

X 

Tv 

rH 

rH 

rH 

CN 

ro 

CO 

tf 

X 

r** 

4 

p* 

o 

in 

CN 

CN 

jo 

in 

4 

r- 

CN 

X 

o 

X 

a 

<r 

X 

X 

rO 

CN 

ro 

CN 

o 

Tv 

CN 

*4* 

1^ 

O' 

o 

X 

in 

X 

r* 

TV 

no 

r— 

n 

CO 

4* 

CN 

CN 

no 

rO 

4 

m 

x. 

X 

H 

ift 

CN 

H 

X 

X 

X 

a> 

■X 

oo 

CN 

O 

ro 

X 

JO 

H 

X 

rH 

«H« 

JO 

X 

JO 

ro 

X 

o 

no 

O 

<r 

'T 

rH 

x> 

■X) 

o 

o 

o 

o 

X 

o 

c 

o 

H 

▼H 

CN 

cO 

X 

CT* 

*o 

r- 

X 

ro 

CN 

X 

r- 

r- 

OD 

CN 

CO 

X 

JO 

X 

ro 

CN 

X 

X 

JO 

X> 

•r 

•* 

N 

£• 

CN 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

H 

rH 

CN 

ro  X 

r- 

O' 

rN 

*0 

rH 

OP 

X 

o 

O' 

JO 

X 

CN 

-t 

n 

r- 

CN 

4 

r»». 

X 

r. 

rH 

rH 

CN 

CN 

ro 

*r 

JO 

r- 

Tv 

JO 

TV 

•H 

0s 

Tv 

;N 

Jv 

JV 

rH 

rH 

rH 

<N 

ro 

CO 

•5- 

X 

T 

x> 

o 

jo 

CO 

no 

p* 

in 

*n 

CN 

JO 

CO 

ro 

CN 

X 

tH 

"0 

X 

o 

ro 

X 

ro 

o 

4* 

4* 

O 

JV 

-o 

Tv 

CN 

X 

■t 

-o 

rN 

CM 

4 

JO 

X) 

n 

p- 

0D 

O 

no 

r- 

rO 

r-i 

ro 

O 

X 

X 

-r 

r- 

n 

r- 

o 

CN 

rH 

■n 

O' 

X 

O 

rH 

X 

rO 

rH 

rf 

cr 

CN 

C 

4“ 

H 

X 

nO 

o 

o 

o 

o 

o 

o 

*H 

«h 

•H 

CN 

r 

X 

•x 

rH 

CN 

CN 

r- 

D 

rH 

rO 

CO 

-M 

X 

T. 

CN 

CN 

X 

o 

Tv 

.n 

n 

•H 

CN 

•H 

rH 

4 

o 

Tv 

0. 

0. 

o 

o 

o 

o 

o 

O 

o 

O 

o 

O 

o 

o 

rH 

H 

CN 

ro 

J0 

r* 

o 

ro 

r* 

CN 

X 

X 

r- 

X 

X 

CN 

JO 

X 

X 

-o 

JO 

X 

rH 

£• 

rH 

H 

rH 

CN 

IN 

ro 

-4* 

X 

T- 

CN 

J0 

Tv 

4 

rH 

TV 

Tv 

N 

JV 

Jv 

rH 

H 

tH 

CN 

ro 

no 

4 

c 

Tv 

X' 

cv 

H 

CN 

00 

CO 

JO 

t> 

no 

rH 

H 

CN 

r*» 

Ov 

X 

C7v 

JO 

Tv 

Tv 

H 

O 

*4* 

rH 

ro 

no 

o 

T> 

X 

O' 

JO 

CN 

Tv 

X 

X 

zr 

sD 

x» 

TV 

o 

hi 

CN 

<4 

r- 

rH 

X 

4 

*T* 

T> 

a> 

X 

JO 

"O 

.*o 

o 

CN 

X 

X 

r- 

tH 

CN 

X 

(N 

r- 

TV 

CN 

rH 

PO 

T\ 

CN 

O 

X 

•4 

*■  i 

X 

X 

rH 

c 

x 

•-I 

H 

*H 

H 

«-i 

CN 

CN 

co 

rr 

00 

o 

3s 

X 

in 

CN 

rH 

X 

X 

CTv 

O 

X 

Jv 

X 

rH 

X 

CN 

CN 

X 

T. 

X 

X 

£ 

X 

*£ 

o 

o 

o 

o 

o 

o 

o 

o 

X 

o 

z> 

o 

o 

rH 

rH 

rH 

CN 

ro 

X 

X 

co 

CO 

Tv 

r- 

o 

ro 

X 

X 

ro 

JO 

X 

vC 

CN 

P- 

* 

tH 

rH 

rH 

CN 

CN 

co 

4* 

X 

r-' 

Tv 

N 

n 

Tv 

4 

rH 

Tv 

•Tv 

fN 

JV 

Tv 

w 

rH 

rH 

rH 

CN 

CO 

ro 

4 

X 

T> 

c 

H 

rH 

<4 

fl 

ro 

-o 

CN 

4 

4 

Tv 

CT> 

•-N 

JO 

rH 

rH 

ro 

CN 

n 

•O 

rs 

n 

ro 

JO 

X 

rH 

c 

4< 

X 

0> 

CN 

CN' 

-4 

CN 

~c 

JO 

no 

•H 

r-» 

rH 

.1) 

00 

£> 

r- 

X 

o 

N 

n 

Tv 

4 

X 

co 

O 

.n 

NJ 

'v-r 

X 

'Tv 

JO 

r- 

X 

r- 

n 

rH 

X 

X 

Tv 

no 

CN 

n 

rH 

■H 

CN 

X 

-** 

4 

rH 

CCl 

«— • 

tH 

CN 

CN 

CN 

CN 

'o 

JO 

X 

X 

O 

CO 

r- 

ro 

o 

T' 

rH 

r- 

CN 

JO 

X 

rH 

no 

X 

n 

X 

^4 

nO 

Ov 

Tv 

X 

D 

JO 

X 

"O 

P- 

4 

c 

o 

t«r 

o 

x 

o 

c 

o 

O 

o 

o 

o 

rH 

•H 

CN 

ro 

n 

in 

X 

X 

rH 

00 

*o 

X 

X* 

X 

H 

n- 

X 

X 

Tv 

4 

-c 

cr 

Tv 

N 

X 

rH 

•H 

•H 

.N 

ro 

co 

4* 

X 

r* 

TV 

CN 

X 

Tv 

t* 

tH 

Tv 

T 

N 

4 

T- 

tH 

t— 1 

rH 

CN 

ro 

'O 

4 

X 

P 

Tv 

4 

TV 

p* 

.H 

CN 

•o 

in 

ro 

CO 

x 

Mr 

r*». 

CO 

CN 

JO 

<r 

<o 

T> 

jc> 

H 

vT 

r- 

•o 

CN 

O' 

*4* 

rH 

Tv 

X 

o 

tH 

n 

n 

*N 

o 

Tv 

c 

T> 

4“ 

o 

t— 1 

■o 

O 

Tv 

CO 

JO 

n 

4 

X 

H 

H 

X 

*o 

X 

ro 

O 

n 

x 

•H 

•4- 

CN 

JO 

rf , 

-n 

O' 

CN 

Tv 

r- 

r- 

o 

C 

O 

X 

n 

\ 

Jv 

£ 

no 

ro 

ro 

ro 

ro 

4* 

4- 

JO 

X 

r- 

cr- 

l 

ro 

H 

in 

X 

r* 

ro 

ro 

Tv 

rH 

CV 

c0 

X 

O 

ro 

CN 

•c 

H 

o 

r- 

m, 

4 

;N 

no 

4 

r— 

O 

o 

x 

x 

o 

o 

3 

X 

D 

3 

o 

rH 

rH 

rH 

CN 

rO 

JO 

av 

rH 

JO 

Tv 

X 

Tv 

Tv 

CN 

X 

T 

J0 

c 

X 

X 

C 

ro 

T 

H 

rH 

-H 

N 

'O 

4* 

X 

P-- 

TV 

<N 

•  o 

Tv 

JO 

rH 

jv 

o 

T» 

T. 

— H 

t— 4 

t-H 

CN 

no 

CO 

JO 

X 

TV 

X> 

vC 

CN 

4* 

n 

00 

in 

iH 

W 

CN 

CO 

X 

X 

X 

X 

o 

ro 

X 

JO 

X 

X 

ro 

X 

JO 

X 

o 

X 

rH 

^O 

iH 

O 

Tv 

T 

rH 

rs 

t> 

CN 

jo 

7v 

4* 

o 

H 

n 

CO 

X 

X 

r+* 

ro 

X 

JO 

ro 

X 

i/o 

.H 

00 

JO 

in 

*o 

ro 

*H 

*4 

rH 

— «* 

Tv 

no 

T> 

ro 

rH 

X 

n 

O 

o 

X) 

in 

■x 

*0 

£> 

r- 

•x 

O' 

x 

ro 

JO 

X 

CN 

r- 

CO 

t— n 

H 

o 

O 

X 

T* 

Tv 

•  1 

X 

X 

tH 

o 

Tv 

4« 

X 

tH 

CN 

-* 

— 

T 

r\ 

X 

o 

c 

O 

o 

O 

o 

o 

x 

H 

H 

«-H 

»-» 

rH 

CN 

CN 

ro 

JO 

X 

X 

CN 

JO 

'Jv 

•JO 

nr 

T> 

ro 

Tv 

o 

X 

-C 

r. 

»H 

X 

X- 

rn 

r 

rH 

rH 

tH 

H 

CN 

ro 

"O 

in 

X 

r* 

o 

'X- 

n 

T> 

JO 

rH 

> 

•X 

✓ 

c* 

rH 

rH 

rH 

CN 

'O 

ro 

n 

o 

p 

o 

•H 

o 

r 

X 

o 

H 

0 

X 

«H 

H 

4 

o 

O' 

JO 

n 

CN 

CN 

CN 

o 

CO 

o 

CJv 

Tv 

ro 

Tv 

^O 

tH 

Tv 

X 

"O 

O 

nO 

T*- 

C 

T> 

X 

rn. 

r^* 

X 

o 

h< 

x 

CN 

JO 

o 

n 

H 

Tv 

x 

*0 

O' 

o 

•.n 

•x 

o> 

H 

X 

X 

rH 

rH 

■T 

X- 

JO 

«H 

'N 

rH 

rH 

Tv 

o 

Tv 

X 

N 

X 

N 

X 

X 

N 

o 

rH 

H 

CN 

CN 

ro 

CO 

JO 

X 

r* 

o 

CN 

JO 

j'. 

LO 

H 

J' 

o 

'O 

JN 

o 

JO 

X 

TV 

•H 

X 

X 

rH 

X 

o 

X 

*4 

X 

rH 

CN 

•H 

H 

O 

X 

rt 

H 

rH 

»h 

H 

H 

<-H 

CN 

CN 

vN  CN 

rO 

X 

X 

rH 

ro 

X 

O 

X 

CN 

o 

rH 

r 

H 

rH 

X 

cr. 

c 

CN 

C 

T 

.N 

tH 

JO 

rH 

H 

•H 

rH 

CN 

CN 

rO 

■41 

JO 

X 

CO 

X 

CN 

X. 

X 

JO 

<H 

Tv 

o 

•“O 

T 

X 

rH 

rH 

rH 

;N 

CN 

^0 

CO 

JO 

vl' 

> 

X 

H 

2 

rH 

TV 

30 

TV 

X 

«r  m 

m 

4 

4 

-H 

O 

> 

X 

X 

T> 

JO 

o 

r 

rH 

fNJ 

rH 

Tv 

JO 

ro 

•T* 

n 

O 

X 

o 

r0 

o 

X 

O' 

r- 

JO 

tH 

T 

X 

0 

o 

< 

rH 

H 

£> 

CN 

r* 

CN 

rH 

rH 

CO 

o 

X 

X 

X 

P** 

CN 

r* 

'O 

X 

C 

J' 

X 

•o 

CN 

0 

i  n 

JO 

X 

<4 

Jv 

X 

'N 

T 

Tv 

rH 

n 

T 

n 

,N 

o 

r 

Ov 

D 

ro 

ro 

tv 

T> 

in 

X 

CO 

r* 

o 

'N 

o 

ro 

o 

X 

0 

H 

ro 

vTV 

•D 

o 

TV 

JO 

CN 

TV 

p* 

X 

-o 

CO 

CN 

r-t 

r- 

H 

H 

a. 

O 

X 

•o 

*o 

m 

tH 

X 

H 

n 

o 

V 

X 

ro 

X 

o 

o 

rO 

rH 

•Tv 

r* 

*n 

~r 

CO 

-o 

CN 

rH 

-H 

rH 

c 

O 

X 

o 

o 

o 

'X 

c 

X 

ON 

.'JO 

X 

X 

P* 

X 

X 

J0 

r> 

'J' 

ro 

'O 

CN 

,N 

>J 

•H 

o 

o 

Cv 

X 

O 

o 

X 

X 

X 

o 

X 

o 

X 

X 

o 

X 

o 

X 

<^« 

o 

o 

o 

o 

o 

o 

o 

x 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

X 

.Hs 

X 

•t 

© 

o 

o 

o 

o 

D 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

X 

o 

X 

o 

o 

X 

X 

y 

< 

o 

o 

H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

r\ 

o 

o 

o 

o 

o 

© 

7* 

JO 

p* 

D 

n 

't 

CO 

CN 

rH 

o 

H 

CN 

rO 

JO 

X 

r- 

it 

JV 

o 

rH 

CN 

■o 

<r 

n 

X 

r^* 

X 

O 

o 

tH 

CN 

'O 

4 

JO 

0 

-J 

X 

i?> 

o 

o 

1 

1 

« 

1 

i 

1 

1 

i 

1 

rH 

•H 

H 

rH 

rH 

H 

rH 

rH 

rH 

CN 

CN 

CN 

N 

CN 

CN 

CN 

>N 

rN 

;N 

ro 

From  the  moments  of  w,  we  determine  the  moments  of  v.  Since 


v  is  the  sum  of  N  of  the  w,  all  of  which  are  i.i.d.,  this  simply 
involves  multiplying  by  N.  Values  of  the  mean,  variance  and  standard 
deviation  of  Wq  (labelled  "single")  and  of  (labelled  "N  Times")  are 
shown  in  Table  5-2. 

Our  initial  approximation  to  the  SNR  required  by  the  GOF  to 
achieve  PD  =  .5  is  given  by  Equation  (3-3).  Substituting  into  this,  we 
have 

N  e[wq(N)]  +  X  >/N.var[wa  (N)]  *  N  Efw^N)]  =  N{e[W()(N-1)]  +  B(a,N)J 

where  the  last  equality  results  from  Equation  (5-13).  Rearranging, 

N  B (a,N)  =  n{e[wq(N)]  -  e[W()(N-1)]  }  +X  >/N.var[w0(N)] 

Using  Equation  (5-11) ,  this  reduces  to 

.(,,II)4*!£S  (5-16) 

VN 

This  is  the  equation  we  solve  to  determine  Y. 

As  before,  we  set  X =4.26489  for  a  false-alarm  probability  of 
10  Using  Tables  5-1  and  5-2  (or  more  finely  granulated  versions  thereof) 
we  then  readily  determine  Y as  a  function  of  N.  This  is  plotted  as 
the  solid  line  in  Figure  5-1. 

Asymptotic  performance  is  not  as  readily  determined  as  for 
the  energy  detector.  However,  we  do  note  that  for  false-alarm  probabilities 
less  than  .084, 

i!T.  Y(N>  > 1  (5*17) 
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To  see  this,  re-write  Equation  (5-16)  as 


VST  B (a,N)  -  -L+  X  yvar[w0(N)]  (5-18) 

yN 

In  the  infinite  limit,  the  right  side  of  this  equation  approaches  Xtt . 
For  Y-l,  we  have  a=^..  For  this  value  of  a,  special  relationships 
involving  the  Gamma  functions  in  Equation  (5-14)  exist.  In  particular, 
from  Equation  (6.1.46)  of  Reference  3, 


lira  1/2  T(N) 
N->»  T(N  4-1/2 


t  We  may  thus  write 

VT  B  (1/2  ,N)  =  T(l/2)  =  y/7 

N->°° 

For  X  =  /6 /it ,  the  solution  to  Equation  (5-18)  as  N+°°  is  thus  a=l/2. 

This  value  of  X  is  equivalent  to  a  of  .084.  For  larger  X  (smaller 

and  more  realistic  PFA) »  a  value  of  u  smaller  than  1/2  (y  greater  than 
unity)  will  be  required  to  satisfy  (5-18),  establishing  our  assertion. 

As  before,  we  now  correct  for  the  error  introduced  in  the 
setting  of  the  threshold.  To  do  this,  we  determine  the  exact  value  of 
A  at  which  the  distribution  of  Vq(N)  is  equal  to  1  -PpA.  Bor  N  =  l,  the 
density  of  v  is  that  of  w,  and  A  is  readily-  found  to  be 

A(N  = 1)  =  -In  PpA  (5-19) 

For  larger  N,  the  density  of  vqOO  is  the  N-fold  convolution  of  the 
density  of  Wq(N).  After  convolving,  this  is  th 'n  integrated  to  obtain 
the  distribution  function  of  Vq(N). 

The  necessary  calculations  are  most  readily  performed  on  a 
computer.  The  density  of  w^(N)  is  generated  by  sampling  over  a  sufficiently 
large  interval.  For  convenience,  the  mean  value  is  removed,  centering 
the  density  at  zero.  An  FFT  then  transforms  the  density  into  the 
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characteristic  function.  Raising  each  term  to  the  N^  power  effects 
the  convolution.  An  inverse  FFT  then  results  in  the  desired  density. 
Numerical  integration  yields  the  distribution,  whence  the  exact  threshold 
(for  given  P^)  can  be  determined. 

From  Equation  (5-12),  the  variance,  and  hence  standard 
deviation,  of  v^(N)  is  known.  By  normalizing  v  by  this  value,  the 
distribution  function  just  calculated  is  expressed  in  terms  of  a  zero 
mean,  unit  variance  random  variable.  By  plotting  this  function  on 
probability  paper,  one  can  determine  the  validity  of  the  normal  form 
used  as  an  approximation.  In  Figure  5-2  is  shown  a  true  normal  distri¬ 
bution  (dashed  line)  and  actual  distributions  of  Vq(N)  for  N  =  16  and 

N=64.  As  expected,  in  the  tails,  the  error  is  noticeable.  At  .9999 

-4 

(P  .  =10  ),  for  N=64,  the  normal  approximation  is  in  error  by  roughly 

FA 

one-third  of  a  standard  deviation. 

Having  determined  the  exact  threshold,  A(N),  we  equate  to  this 
the  mean  of  v^(N).  Solving  for  y,  we  obtain  the  first  correction  to 
our  initial  approximation.  This  is  shown  as  the  dashed  curve  in  Figure  5-1. 

The  second  source  of  error  arises  from  equating  the  mean  to 
the  median  of  v^(N) .  Using  the  same  approach  as  for  v^,  we  may  find  the 
density  of  v^(N)  as  the  N-fold  convolution  of  the  density  of  w^(N). 

The  procedure  is,  however,  not  as  straightforward  since  the  density  of 
Wp(N)  is  a  function  of  a(or,  equivalently,  of  Y).  A  value  of  Y  is 
chosen,  the  density  of  v^(N)  is  then  calculated  and  integrated.  The 
fiftieth  percentile  point  is  compared  to  the  exact  threshold  A(N) .  Y 
is  varied  until  the  two  coincide.  This  value  of  signal-to-noise  is 


the  exact  solution. 
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For  N  = 1 ,  the  exact  value  of  Y  is  already  known  to  be  11.9  dB. 
For  N  =  2  and  N  =  8,  exact  values  (of  Y=9.1  dB  and  Y=5.7  dB)  were  found. 

As  the  corrections  rapidly  become  small,  further  values  were  not  found. 

The  dashed  and  dotted  curve  in  Figure  5-1,  fitted  through  these  points 
to  the  asymptotic  result,  thus  represents  the  performance  of  the  "greatest 


of"  receiver. 
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SECTION  6 

COMPARISONS  AND  CCNCLUSIONS 

As  explained  in  Section  2.4,  we  are  comparing  an  interceptor 
operating  at  a  per-decision  false-alarm  probability  of  10  against  the 
intended  receiver  operating  at  10  ^ .  The  performance  of  the  intended 
receiver  operating  at  this  is  shown  in  Figure  4-4  of  Reference  1. 

From  this,  we  extract  the  values  of  SNR  required  by  the  intended  receiver 
to  have  a  fifty  percent  detection  probability  as  a  function  of  N.  These 
are  shown  in  Figure  6-1.  Our  comparison  consists  of  simply  subtracting 
the  required  SNR  for  the  ED  and  GOF  from  that  for  the  intended  receiver. 
This  yields  the  desired  recognition  differential,  which  is  shown  in 
Figure  6-2. 

The  negative  values  of  recognition  differential  observed  for 
small  N  reflect  the  difference  in  false  alarm  levels.  If  both  the 
intended  receiver  and  the  interceptors  were  operating  at  the  same  false 
alarm  level,  then  the  recognition  differential  would  be  zero  at  N  = 1 
and  strictly  positive  thereafter. 

The  ED  is  viewed  as  representing  an  upper  bound  on  recognition 
differential,  similarly,  the  GOF  is  a  lower  bound.  Against  any  other 
form  of  intercept  receiver,  the  recognition  differential  inherent  in  the 
chosen  message  structure  will  be  bounded  by  these  two  curves.  As  the 
difference  between  the  two  is  not  very  great,  we  conclude  that  there 
is  no  particular  need  to  analyze  in  detail  any  of  the  other  possible 
strategies  available  to  an  interceptor. 
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FIGURE  6-2  RECOGNITION  DIFFERENTIAL 
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APPENDIX  A 


SOME  SUMMATIONS  INVOLVING  BINOMIAL  COEFFICIENTS 


In  this  appendix  we  present  the  details  behind  the  derivation 

of  the  moments  of  Wq(N)  and  w^(N)  for  the  "greatest  of"  receiver.  First, 

2 

we  determine  the  expected  value  of  w  and  w  under  each  hypothesis.  We  . 
then  simplify  our  results  and  determine  expressions  for  the  mean  and 
variance.  By  definition. 


eCw0W3  =  /  wf0  N(w)dw 

0 

(A-l) 

E[w02(N)]  =/  w2  fQ  N(w)  dw 

0 

From  Equation  (5-8) ,  we  have 

(A-2) 

W”>  - 

(A-3) 

which  we  re-write  as 

f»  „(»)  -  »  ^  Ck) 

k=0 


by  use  of  a  binomial  expansion.  Utilizing  a  basic  property  of  binomial 
coefficients,  viz. 


(A-4) 


we  have 


N-l  /  N  \  ,  k+i 

fo,,«  -E  U 

k=0 


which,  via  a  change  in  index,  becomes 


Substituting  this  into  Equation  (A-l)  and  interchanging  the  summation 
and  integration,  we  have 


N  aN  \  « 

■[wo0O]  “53  'k'  (_1)k"llc  J  we_kWdw 

k=l  0 

2 

The  integral  is  trivially  shown  to  evaluate  to  1/k  ,  hence, 

N  /N, 


EC“0<N)3  (k)  (-»k'1(f) 


(A -6) 


k=0 


2 

The  evaluation  of  e[wq  (N)]  is  almost  identical.  The  integral 
in  this  case  is 


2  -kw,  2 

-  e  dw  ■*»  — r 

k 


so  that 


/"2 

e["o2<n)1  -  0  <-1>k‘1  (i) 


(A-7) 


with 

and 

Then, 


k=l 


We  now  consider  H^.  Using  Equation  (5-10),  we  write 

fl,N^  "  "A  46 

11  (H-l)  (l-e"w)E_2  e”(° +  Dw 
B1  ae-^Cl  -  e"”)"'1 


00 

Efwj^CN)]  -  e[w0(N-1)]  +  I  w  (B-A)dw 

0  00 

e£w^(n)]  =  E[wo^(N-l)]  +  J  w^  (B-A)dw 

0 

We  expand  A  and  B . 

_  /N“l\  /  ,\k,  -Wvk  -avs 

B  «  a  >  ,  (  J  (-1)  (e  )  e 


(A-8) 

(A-9) 


k=0 


A-ff-l)g  (N_2)  (-l)k(e-)k+1  e' 


aw 


k“0 


A-2 


The  expression  for  A  is  re-written  employing  the  same  technique  used  to 
change  Equation  (A-3)  into  (A -5) . 

A“Z  O  s-0” 

k=l 

We  note  that  for  k  =  0,  the  summand  of  A  is  zero,  hence  we  may  re-define 
the  summation  for  A  as  commencing  at  k  =  0.  This  allows  us  to  combine 
A  and  B  into  a  single  summation: 

B-A-jrf  f  '1)  (-l)k  [(k  +  a)e-(k  +  “>”]  (A-10) 

k=0  k 

Substitution  of  Equation  (A-10)  into  (A-8)  and  (A-9)  followed  by  an  inter¬ 
change  of  the  summation  and  integration  yields  results  very  similar  in 
form  to  those  stated  in  Equations  (A-6)  and  (A-7) . 


N-l  N. 

m  1  y  |  /  1  \ 

£[^00]  =E[w0(N-l)]  +2 3  ( 

k=0  '  1 

t )  <-1>k  (£)  . 

(A-ll) 

N-l 

B.l.  ,  .2 

E[Wl2(N)]  =E[wo2(N-1)]  +2^  ( 

k=0 

:vi  ->*  <*> 

(A-12) 

Note  that  the  summations  now  start 

at  k  =  0,  not  k  =  1. 

Thus,  as  a-*0. 

the  leading  term  of  the  sums  approaches  ® .  This  is  to  be  expected  as 


«  =  0  corresponds  to  an  infinite  signal -to -noise  ratio,  at  which  the 
expected  value  of  the  "greatest  of"  should  be  infinite. 


We  now  simplify  the  summations  for  the  mean  and  variance  of 


Wq(N)  and  the  mean  of  w^N). 

The  mean  of  Wq(N)  Is  defined  by  Equation  (A-6).  We  use  an 
inductive  approach  to  express  it  more  compactly.  Let  f  represent 
the  summation  up  to  n,  i.e.. 


(A-13) 


By  suitable  manipulation,  this  is  expressed  in  terms  of  f  whence  the 
behavior  as  a  function  of  n  is  observed.  We  use  the  readily  proven 
recurrence  relationship, 


which  allows  us  to  write 


f 


n 


(A-14) 


In  the  first  summation,  the  k«n  term  is  zero  since,  by  definition. 


C)  ■  0 


Thus,  the  first  summation  may  be  written  as  from  k  =  l  to  n  -1, 
is  directly  seen  to  be  f  Using  Equation  (A -4) ,  the  second 
re-written  as 


which 
sum  is 


I 

n 


2  <-l)k'X 


k=l 


We  then  have 


> 


f 

n 


fn-l+nE  <*l> 
k-1 


k-1 


A-4 


However,  since 


[l  +  M)]n-£  <-»k  (")  .0 


k=0 

we  have  the  recurrence  relationship 


f  =  f  .  +  - 
n  n-1  n 


Noting  that  fj=l,  we  finally  obtain  the  closed  form  expression 


n 


(A-15) 


k=l 


A  similar  procedure  is  now  applied  to  the  expected  value  of 


Wq  (N) .  Defining 


8--S(-l)k'1  Oft) 


we  readily  establish  that 


As  g^  =  0  and  g^  =  1 ,  we  have 


3  8n-l 

+  -  f 
n  n 

have 

n  _ 

n  k 

E  k 

*k 

« 

M 

M 

_1 

m 

k=l 

k=l  m=l 

of  V 

(N) 

is  directly 

expressible 

(A -16) 


(A-17) 


and 


V 


var  [w0(N)]  -  2gR  -  fj 


(A-18) 


Expanding  each  of  these  as  a  double  summation, 

k«=l  m=l  k=l  m^l 
H  r  k  .  N  -| 

-E  ■  E  il 

k*l  **m=l  m**l  J 


A -5 


This  is  equivalent  to 


where 


m  <  k 


m  >  k 


Clearly,  all  the  terms  cancel  except  for  those  for  which  k=m.  Thus, 


var  [w0(N)] 


(A-19) 


Finally,  we  examine  the  mean  of  w^(N).  We  consider  the  summation  in 
Equation  (A-ll)  and  claim  that  it  is  equivalent  to  a  Beta  function. 
This  function  is  defined  by 


B(N 

Expanding  the  integrand. 


, a)  «=  B  (a,N)  A  J  a-t)N“1t“”1  dt 


B(N,a) 


and  interchanging  the  order  of  operation, 


■s 


or  B(N,a)  =  2  (-l)k  ^  X)  kfcl 

k=0  '  K  ' 

proving  our  claim.  We  may  thus  write 

E[Wi<N)]  =  E[w0(N-l)j  +  B(N,o) 


(A-20) 


A-6 


(A-21) 


APPENDIX  F 


PROCESSOR  SIZING 


The  processing  required  to  implement  the  double-level  soft 
decision  decoder  technique  proposed  for  the  HIDAR  Receiver  has  been  analyzed 
in  order  to  provide  a  measure  of  assurance  that  the  approach  could  be 
realized  in  hardware.  The  candidate  processor  selected  is  the  Data 
Dependent  Signal  Processor  (DDSP)  recently  developed  by  Steir.  Associates. 
This  processor  is  soon  to  be  announced  as  model  AR-10/10  of  the  AR-10 
family  of  modular  high-speed  digital  processors.  While  the  DDSP  has 
certain  special  features  that  make  the  proposed  processing  more  efficient 
and  relatively  easy  to  program,  it  is  not  unfair  to  say  that  the  DDSP 
represents  the  general  level  of  processing  ability  that  is  available  with 
a  typical  modern  signal  processing  computer. 

Figure  F-l  is  a  flowchart  illustrating  the  main  processing  tasks 
that  need  to  be  performed  to  implement  the  double-level  decoding.  The 
flowchart  is  somewhat  simplified,  neglecting  the  facts  that  1)  double¬ 
level  accumulation  cannot  begin  until  the  start  of  message  condition  is 
detected,  2)  first  pass  is  simpler  (there  are  no  previous  accumulations), 
and  3)  the  last  pass  is  simpler  since  only  one  value  of  K  (see  below)  is 
required.  For  timing  analysis  it  is  enough  to  recognize  that  as  these 
special  situations  are  all  simpler,  the  maximum  burdens  occur  during  a 
typical  pass  as  the  message  is  being  received. 

The  first  task,  accepting  the  input  samples,  may  be  made  nearly 
transparent  if  special  hardware  is  provided  to  put  the  incoming  data  in 
place  in  the  working  memories  of  the  DDSP  using  a  background  DMA  data 


(U)  BASIC  DECODER  FLOW  (U) 
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transfer.  Since  this  technique  is  unique  to  the  DDSP,  we  will  take  a 
worse-case  assumption  that  such  hardware  is  not  present  requiring  the 
program  to  sense  that  the  data  is  ready  on  a  sample-by-sample  basis. 

Once  the  processing  of  a  block  of  data  is  started,  it  will  be  necessary 
to  interrupt  processing  at  regular  intervals  to  save  input  data  and  store 
it  away  in  a  software  FIFO  (first-in,  first-out)  memory  for  later  use. 

When  the  processing  for  the  previous  block  is  complete,  the  program  must 
first  empty  the  accumulated  samples  in  the  FIFO  before  sensing  new  samples 
at  the  input.  Such  a  routine,  coded  for  the  DDSP  as  part  of  a  previously 
developed  program,  times  out  at  approximately  1.5  us  per  sample.  Since 
there  are  assumed  to  be  fo’-ty  tones  per  subband,  the  number  of  samples 
per  pass  is  40x7=280  requiring  280  x  1.5  us  s  .4  ms. 

Once  the  data  is  taken,  »  complex  FFT  must  be  performed  using 
a  number  of  points  equal  to  the  incoming  sample  set  rounded  up  to  the  next 
binary  increment.  The  FFT  has  been  coded  for  the  DDSP,  requiring  12 
instructions  (1.62  ms)  for  the  butterfly  calculation  that  is  the  heart 
of  the  process.  With  approximately  10%  additional  processing  to  set  up 
the  parameters  and  indexes  for  each  stage,  the  worst  case  (512  points) 
total  time  required  is  (1 .62x2304)1 ,1 ~4 .1  ms. 

After  the  FFT,  the  program  operates  on  a  subband  basis,  calculating 
2 

the  (magnitude)  for  each  frequency  (using  only  odd  or  even),  correlating 
again  all  possible  codes  and  updating  the  double-level  accumulators  and 
vectors . 

Of  these  processes,  forming  the  square  of  the  magnitude  is  most 
straightforward.  This  process  requires  two  multiplications  (two  instruction 
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times  each)  plus  a  summation  with  one  or  two  more  instructions  to  fetch 
and  restore  the  data.  By  overlaying  the  processing  for  one  frequency  with 
the  next,  the  total  processing  can  be  reduced  to  five  instructions  per  sample 
for  a  total  burden  of  5  instructions  x  20  frequencies  plus  20  instructions 
to  initialize  the  loop.  The  total  burden  is  120  x  135  ns  x  7~.l  ms. 

The  correlation  process  is  somewhat  more  complex  since  it 
consists  of  selected  accumulations  of  ten  of  the  twenty  input  magnitudes 
where  32  different  patterns  for  selection  are  used.  The  innermost  part 
of  this  program  consists  of  a  subroutine  that  reads  in  the  current  pattern 
word,  examines  the  next  bit  and,  if  it  is  ONE,  adds  a  corresponding 
magnitude  into  an  accumulator.  This  routine  is  estimated  to  require  4 
instructions  if  the  pattern  bit  is  ZERO  or  6  instructions  if  the  pattern 
bit  is  ONE.  Since  half  the  bits  are  ONE  and  half  ZERO,  the  average  is 
five  instructions.  With  twenty  pattern  bits  per  code,  32  codes,  five 
instructions  per  code,  135  ns  per  instruction  and  7  subbands,  the  total, 
allotting  10  percent  for  overhead,  is  (20  x  32  x  5  x  135ns  x  7)1.1 ^3.3  ms. 

Most  complex  of  all  is  the  double-level  decoding.  The  innermost 
loop  of  this  process  requires  the  addition  of  the  ith  accumulator  with 
the  jth  correlation  sura  where  the  sum  of  i  and  j  is  K(modulo  32).  After 
the  32  summations  that  meet  this  criteria  are  formed  the  largest  is  selected 
and  the  new  sum  replaces  the  old  and  the  number  j  is  added  to  the  decision 
vector.  This  process  is  repeated  for  all  32  possible  values  of  K. 

The  program  to  perform  this  inner  loop  must  first  calculate 
the  j  associated  with  each  i,  fetch  the  corresponding  accumulator  and 
correlation  sum,  add  them,  test  to  see  if  it  is  larger  than  previous 
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sums  and,  if  so,  save  the  sum,  i  and  j.  This  routine  has  been  tentatively 
programmed,  requiring  12  instructions.  Since  it  is  in  a  double  loop 
each  repeated  32  times,  the  total  requirement  is  {  (12  instructions  x  32)1.1 
x  32  x  135  ns /instruct ion  x  7  subbands}  1.1~14.1  ms. 

Since  the  total  burden  is  the  sum  of  the  individual  requirements 
it  is  estimated  that  the  total  processing  is  approximately  .4ms  +4.1  ms 
+  .1  ms  +  3.3  ms  +  14.1  ms  =  22.0  ms.  This  processing  is  repeated  every 
30  ms.  Thus,  the  burden  may  also  be  expressed  as  22/30  ^  73%  of  available 
processing  time,  a  figure  that  provides  an  acceptable  level  of  margin 
allowing  for  the  growth  and  additional  tasks  that  usually  must  be  added 
as  a  system  moves  from  the  analysis  stages  to  detailed  implementation. 


